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ABSTRACT 


This report summarizes work to date performed under NASA Grant 
NGR 34-002-195, entitled A Theoretical Analysis of the Current -Voltage 
Characteristics of Solar Cells. The objective of the research is to 
provide a theoretical understanding of the operation of solar cells 
so they can be designed for maximum efficiency. These objectives are 
being pursued through the application of a general computer device 
analysis program to the study of solar cell properties. During the 
past twelve months considerable progress has been made toward these 
objectives. Specifically an analysis program for solar cells has been 
developed which provides for non-uniform doping profiles , generation- 
recombination effects, field and doping dependent mobility, and optical 
generation. Preliminary results wirh the analysis program indicates 
that the predicted efficiencies agree reasonably well with experimental 
efficiencies for silicon solar cells. 

This report summarizes the major work to date on the current- 
voltage characteristics and efficiencies of solar cells. For one solar 
cell structure detailed curves are presented which include carrier 
densities, current densities, potential, and quasi-Fermi levels at 
different voltage levels both with and without optically generated 
carriers (AMO conditions) . In addition some results are presented 
concerning the influence of various parameter variations such as 
lifetime, cell thickness, and high-low junction width on solar cell 
performance . ' 



I. COMPUTER ANALYSIS PROGRAM 


A considerable part of the effort during this first year's 
research has been devoted to adapting a computer device analysis 
program to the analysis of solar cells. The present computer 
analysis program is able to account for the following physical effects 
in a semiconductor device: 

1. Drift and diffusion currents 

2. Non-uniform doping profiles 

3. Bulk generation-recombination effects 

4. Field dependent and doping dependent mobility 

5. Optical carrier generation 

6 . Surface recombination 

The effects above which have been included in the analysis program 
during the present research work are non-uniform doping profiles s 
optical carrier generation , and surface recombination. 

Two papers have been prepared for submission for publication. 

One of these discusses, the computer analysis program and the mathematical 
techniques used to solve the nonlinear semiconductor device equations , 
while the other discusses the studies which have been carried out on 
optical carrier generation and the effects of ant ire fleeting layers on 
optical carrier generation in silicon. Copies of these manuscripts are 
included in Appendices A and B. These provide detailed summaries of 
the work performed on developing the computer analysis programs. 
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II. DETAILED CHARACTERISTICS AND PROPERTIES OF 
n + -p-p + SOLAR CELLS 

A considerable number of detailed calculations have been made of the 
electrostatic potential} quasi— Fermi levels, carrier densities, and current 
densities in solar cells* Both the conventional n -p cell and the back 
surface field (or high- low junction) n -p-p cell, as shown in cross 
section in Figure 1, have been analyzed. More effort has been directed 

t 4* 

toward the n -p-p structure than the conventional cell because of 
questions concerning the enhanced efficiency of the n -p-p cell. 

Since the device analysis program solves for the internal potential 
and quasi-Fermi levels, it provides a detailed look at the internal 
operation of the solar cell. A large amount of data is generated on each 
device studied. The results presented in this section show detailed 
internal characteristics for one particular solar cell under both conditions 
of no illumination and under the full solar spectrum illumination (AMO ) . 
Material and dimensional properties of the device for which detailed plots 
are presented are given in Table I. The lifetime of 100 ysec in the 

Table I . Material and dimensional properties 


Overall cell thickness 

250 yM 

n + thickness 

0.25 yM 

p + thickness 

0.50 yM 

n surface concentration 

10 2 °/cm 3 

p doping concentration CIO £2* cm) 

1.3xL© 15 

p + doping concentration 

10 18 /om 3 

minority carrier lifetime t 

no 

100 usee 

T 

Ocl ysec 

po 



0.1 ysec 
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p-region corresponds to a diffusion length of about 588 pM or a little 
more than twice the cell thickness. The first sequence of Figures 2-8 
show the device characteristics without illumination. Figure 2 shows 
the electrostatic potential across the device with 2(a) showing the 
entire device length while 2(b) and 2(c) shows expanded views about 
the p-n + and p"*"-p junctions respectively. The various curves correspond 
to applied voltages ranging from 0.1 to 0.7 volts. The expanded views 
are needed to see details about the junctions because of their small 
thickness. Of special interest in this figure is the reduction in 
potential across the high-low junction as seen in Figure 2(c) for 
large values of applied voltage. This reduction begins to occur when 
the minority carrier density exceeds the doping in the p-region, _i 
when. high injection occurs. 

Figures 3(a), 3(b), and 3(c) show complete and expanded plots of 
electron quasi-Fermi level in the device. The expanded plots illustrate 
that the quasi-Fermi level is essentially constant across both junction 
space charge regions as. is assumed in most first order analytical device 
calculations. This first order assumption is not made in this calculation 
but it is seen that relatively constant quasi-Fermi levels in the space 
charge regions do in fact result. The major feature of these plots is 
the almost abrupt change in as seen in Figure 3(c) at the high-low 
junction. This change in slope is required because of the essentially 
constant current density across the high-low junction coupled with the 
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constant in this region as seen from Figure 4(a). Changes in $ 
only occur very near the surface , and these large changes only occur 
because of the value of surface recombination used in the calculations 
(infinite value used). 

Figures 5(a), 5(b), and 5(c) show electron concentration for the 

complete device and expanded views near the two surfaces for various 

applied voltage values. The large changes in carrier density across 

the high-low junction as seen in Figure 5(c) illustrate the minority 

carrier reflecting properties of the high-low junction. The log scales 

used to plot electron density are needed to cover the large range of 

densities encountered as the applied voltage ranges from 0.1 volt to 

0.7 volt; however, these log scales tend to distort somewhat the 

carrier density plots. For example the electron concentrations in the 

p + region as shown in Figure 5(c) from 0 to 0.5 pM vary approximately 

linearly with distance. This linear variation plots on a log scale 

as a very rapidly varying function near aero and a slowly varying function 

as 0.5 pM is approached. The very rapid decrease in electron concentration 

near the back contact is due to the assumption of an ohmic contact. This 

does, not allow the carrier density to increase from its equilibrium value 

2 3 

which is slightly above 10 /cm , The change in electron concentration 
across the high-low junction is related to the high-low junction potential 
as 

n + = expC-qV^/kT) . (1) 

The high-low junction potential is given as 

.N + 

V M = — An r-JpAl + n/K )]. 
hi q N P P 


( 2 ) 
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The combination of Eqs. (1) and (2) gives 

N n 

n • + = n ( 1 + . ( 3 ) 

P P V N p 

This relationship, as discussed in Appendix C, is followed very accurately 

for the high- low junctions analyzed in this work. 

Hole concentration curves are shown in Figures 6(a), 6(b), and 6(c) 

for both the complete device and expanded views about the two junctions. 

The hole density in the center of the device is seen to be constant and 

essentially independent of applied voltage except for the highest 

voltage values which, correspond to high injection in the center p-region. 

Of special interest are the hole densities in Figure 6(c) for the 

n + surface region. This illustrates the injection of holes into the 

n + region (x<0.25 jjM) . The hole density exhibits a rapid decrease near 

x=0 due to the large surface recombination velocity assumed at the surface. 

For low values of surface recombination the hole density will not be forced 

2 3 

to the low values (below 10 /cm ) shown in Figure 6(c) . 

Electron and hole current 'densities as functions of distance at the 
various applied voltage values are shown in Figures 7 and 8. In this case 
expanded views are shown only for the surface n + -p junction since the 
current densities are essentially constant through the p + region. As can 
be seen in Figures 7(b) and 8(b) at low forward voltages there is a large 
change in current density in the vicinity of the n -p junction 
(x < 0.5 jjM ) . This occurs because a large part of the current density 
consists of recombination current from within the junction depletion 
region at low applied voltages. 
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A set of curves similar to those of Figures 2-8 are shown in 
Figures 9-15 for the same solar cell structure but including pair 
generation due to AMO illumination. The assumed surface condition is 
that of a polished surface with no antireflecting layer but reflection 
from the surface is taken into account. For this sequence of figures, 
the terminal voltage ranges from 0.6 volt to 0.0 volt with curves shown 
for 0,1 volt steps. The irradiated surface is at x = 250 pM. The 
open-circuit voltage condition occurs between the 0.5 volt and 0.6 volt 
curves. For all curves except the current density curves, the curves 
do not exhibit any unusual features near the open-circuit voltage 
condition. 

The electrostatic potential curves of Figure 9 are for all practical 
purposes identical to those of Figure 2 in the absence of illumination. 

The potential distribution within the device is determined by the terminal 
voltage and is independent of the optical illumination. 

„The variation,, in the quasi-Fermi level can best be illustrated in 

conjunction with the carrier density curves. The illumination produces 

an excess electron and hole density which is quite evident in Figures 12 and 

13. This excess carrier density is also indicated by the shift of the quasi- 

Fermi potentials. This shift is most evident for electrons in the p-region 

(<£ of Figures 10(a) and 10(b)) and holes in the n region (<j> of 
n p 

Figure 11(b)). At lower bias the electron concentration in the p-region 
(Figure 12(a)) is greatly reduced near the n + -p junction due to the 
collection of the electrons by this junction yielding the short circuit 
current. It can be noticed, however, from Figure 12(b) that the decrease 
in electron concentration near the junction is limited by the saturation 
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drift velocity. The minimum electron density of around 2 x 10 /cm" 

(near x = 0.5 yM in Figure 12(b)) is needed at the saturated carrier 
drift velocity to support the collected photocurrent. This is also 
present in Figure 13(b) which indicates that the saturation velocity 
is reached for holes . This is seen as the fairly flat region of hole 
density below x = 0.5 yM. The surface variation of p^* and <(> , in 
Figures 13(b) and 11(b) also indicate the directional change in J due 

r 

to surface recombination . The collection of the optically generated 
holes is evidenced by the small reverse curve in Figure 13(b) for 
voltages less than 0.4 v. Had velocity saturation not set in this 
would have been a much deeper curve. It is none the less quite evident 
that first order diffusion methods cannot properly model the case of 
high optical currents. 

The effects of the various .-saturation mechanisms point out the 

need for future calculations to include the full generation term instead 

of the present methods of scaling up the results from the case of lower 

generation due to the absence of an antireflection layer. Included in 

these effects is that of high injection due to the carriers. It can be 

seen from Figure 12(a) that high injection (due to optical generation 

alone) is approached in the bulk p-region as the electron density 
13 3 

approaches 10 /cm . For the situation of higher bulk generation rates 
it is quite likely that the high injection condition in the bulk would 
be met. Among, other things this would alter the reflecting characteristics 
of the high-low junction. These conditions of drift velocity saturation 
and high injection could have significant effects on the terminal 
characteristics which would not be evident from calculations based alone 
on the lower generation rates. 
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The current density curves of Figures 14 and 15 are, as before, 
plotted in terms of absolute magnitude . The dips in these curves are in . 
fact directional changes so the current does pass through zero at these 
points. There is no directional change in electron current for high 
bias (V = 0.6 v, Figure 14(a)) due to the overriding effect of the 

cl 

forward p-n junction current* At lower bias however ? the predominant 
electron current consists of collection of optically generated carriers 
by the p-n + junction. The 0.5 v curve in Figure 14(a) represents a 
transitional region around the open circuit condition. For situations 
where the effective surface recombination velocity of the high-low 
junction is much lower, it has been observed that the directional change 
due to carriers supporting back contact recombination occurs much 
closer to the high-low junction than that illustrated in Figure 14(a). 

The hole current density of Figure 15 illustrates a similar 
transitional curve but at a higher bias voltage (0.6 v) due to the lower 
forward pn junction hole current. The current lost due to surface 
recombination is quite evident from Figure 15(b), where a reversal of 
the current direction near the surface indicates current flow to the 
surface due to surface recombination. Realise that for the lower bias 
(less than 0.6 v) the predominant hole current consists of the optical 
component (except at the surface). 

The terminal characteristics of this solar cell are shown in 
Figures 16, 17 and 18. The J-V characteristic in the absence of light 
is shown in Figure 16. This curve is typical of most silicon diodes. 

At voltages ranging from about 0.35 volt, to about 0.55 volt the current 

density varies approximately as exp(qV/kT). At lower and higher voltages 
the current changes less rapidly with voltage approaching the approximate 
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variation of exp(qV/2kT). The deviation from ideal behavior at low 
voltages arises from recombination of carriers within the junction 
depletion region while the high voltage deviation occurs because of 
high injection into the p-region. 

The terminal J-V characteristic under illumination is seen in 
Figure 17. From this the short circuit current is seen to be about 
34 ma and the open circuit voltage about 0.58 volt. The efficiency 
vs voltage characteristic is seen in Figure 18 where the peak efficiency 
is slightly below 12%. This is the calculated efficiency without an 
antireflecting layer. The efficiency of such a cell with an anti- 
reflecting layer would be expected to be about 40% larger or between 
16% and 17%. This calculated efficiency is significantly larger than 
that observed in commercial 10 Q*cm solar cells. There are several 
possible reasons for this. First the lifetime of 100 ysec used in the 
p-region may be larger than that obtainable in present solar cells. 
Lifetime values of 10 ysec have frequently been quoted as typical 
lifetimes in 10 £3* cm solar cells. Calculations for lower lifetime 
cells are given in the next section. The lifetime of 0.1 ysec used 
in the n + region may again be larger than actually obtained in commercail 
cells. However, very small values of lifetime are required in the n + 
layer to significantly reduce the efficiency because of the very thin 
n + region. It was also found that the calculated generation rate used 
in the presented data was about 10% high. As a result of this the 
efficiency results actually lie between that of a bare surface and that 
with a reflecting layer. •• Including this factor would lower the efficiency 
for the .-.cell .with: the antireflecting layer to between 15 and 16%. The 
optical generation parameter presented in Appendix B represent the 
corrected generation rate. 
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There are other effects which may reduce the efficiency in real 
solar cells. The ohmic resistance due to a ’two-dimensional contact 
geometry reduces the curve factor decreasing the efficiency. Also 
heavy doping (degeneracy) effects in the emitter are not accounted 
for. These effects can enhance the loss of carriers at the surface 
and decrease the efficiency. These effects are to be incorporated 
into the theoretical analysis during the continued research work. 
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Figure 2c. Electrostatic potential variation about the high-loy? 
junction Cno illumination) . 
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Figure 4a. Hole Quasi-Fermi potential variation throughout the entire 
device (no illumination). 
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Figure 5a. Electron concentration throughout the entire device 
(no illumination) . 
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Figure 6a. Hole concentration throughout the entire device 
(no illumination) . 
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Figure 6b. Hole concentration about the n + -p junction (no illumination) j 
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Figure 6c. Hole concentration about the high-low junction 
(no illumination). 
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Figure 7b. Electron current density variation about the n + -p junqtion 
(.no illumination). 
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Figure. 9a. Electrostatic potential variation throughout the entire 
device (AMO illumination) . 
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Figure 9c. 


Electrostatic potential variation about the high-low 
junction (AMO' illumination) . 
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Figure 10b. Electron Quasi-Fermi potential variation about the p-n + 
junction (AMO illumination). 
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Figure lib. Hole Quasi-Fermi Potential variation about the p-n + 
junction (AMO illumination). 
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Figure 12a. Electron concentration throughout the entire device 
(AMO illumination). 
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Figure 12c. Electron concentration about the high-low junction 
(AMO illumination).. 
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Figure 13a. Hole concentration throughout the entire device (AMO 
illumination) . 
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Figure 13c. Hole concentration about the high-low junction 
■ (AMO illumination). 
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'igure 14b. Electron current density variation about the' n -p junction 
(AMO illumination). 
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Figure 15a. Hole current density variation throughout the entire 
device (AMO illumination) . 
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Figure 15b. Hole current density variation about the high-low 
junction (AMO illumination). 
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Figure 17 Terminal J-V characteristics of the n + -p-p solar cell under AMO illumination* 








III. EFFECT OF HIGH-LOW JUNCTION ON SOLAR CELL 
PROPERTIES AND EFFICIENCY 
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In addition to developing a solar cell analysis program, one 
of the major efforts during this year's research has been the 
attempt to better understand the role of a p + -p or high-low junction 
at the back surface of a solar cell. In general, the computer 
calculations have shown that the first order models for high-low 
junction behavior are very closely satisfied by the junctions studied 
in this work. A comparison of the numerical results with first order 
high-low junction equations is given in Appendix C. 

The high-low junction has been found to act at low voltages 
primarily as a minority carrier reflecting back contact giving a small 
effective surface recombination velocity at the boundary of the p-p + 
junction. The curves of electron concentration in Figure 5 show the 
large change in electron concentration occurring at the high-low 
junction. 

Several calculations have been made to explore the effect of a 
back surface high-low junction on overall efficiency. Figures 19, 20 
and 21 show one series of modifications of the back p + contact. 

Curves are shown for three devices: (a) ohmic back contact, (b) 0.5 yM 

thick back p region and (c) 5 yM thick back p region. Other parameters 
for the devices are as given in Table I . Figure 19 shows a very 
significant decrease in the current density at a constant voltage for 
the devices with p + back contacts. The decrease in current is due to 
the decreased flow of injected electrons to the back ohmic contact. 

The 5 yM p + device shows only a slight improvement over the 0.5 yM device 



52 


The 5 yM and 0.5 yM devices have back boundary effective surface 
recombination velocities of 1.1 cm/sec and 204 cm/ sec respectively. 

The value for the 0 • 5 yM device is already so low that it acts almost 
as an ideal minority carrier reflecting back contact. 

The current density vs_ voltage characteristics for the three 
devices are shown in Figure 20. The ohmic back contact device is 
seen to have a reduced short circuit current and a reduced open 
circuit voltage . The reduced short circuit current ' arises from the 
loss of optically generated carriers at the back contact. The 
reduced open circuit voltage arises from the combined effects of a 
smaller photocurrent and a larger forward diode current as indicated 
in Figure 19. The efficiency curves of Figure 21 illustrate the 
improved efficiency of the high-low junction solar cell. There is 
little change in peak efficiency between the 0.5 yM and 5 yM cells . 

Calculations similar to those described above were also made for 
250 yM thick devices with 10 ysec minority carrier lifetimes in the 
p-region. These calculations showed a considerably smaller difference 
in peak efficiency between the high-low junction cells and ohmic back 
contact cells. For 10 ysec devices, the efficiency of the high-lpw 
junction device was only about 0.5% above that of the ohmic contact 
device. These results are to be expected on theoretical grounds. 

For a 10 ysec lifetime * the diffusion length in the p-layer is about 
18Q yM which is less, than the 250 yM thickness. On the other hand 
for a lifetime ©f 100 ysec, the diffusion length is about 580 yM or 
slightly more than twice the cell thickness. As the cell thickness 
becomes large compared with a diffusion length, there can be little 
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interaction between the minority carriers and the high-low junction and 
the behavior of the p + device must approach that of the ohmic back 
contact device for small diffusion lengths. 

The voltage change across a high-low junction based upop simple 
theoretical arguments is: related, tO’Wthe -carrier densities from . 

Eq. (2) as 


kT 


N + 

V.. = — in E-g-V Cl + n /N )] 
hi q N p p 

F 


= V. 


kT 


, . in (1 + n /N ) , 

hlo q p p 


(4) 

(5) 


where V hlQ is the equilibrium high -low junction potential. Equation 5 
predicts a reduction in the high-low junction voltage due to the electron 
concentration in the p-region, n^. Limiting cases of Eq. (5) are 


r 


V St 

hi 




v j> 

hlo q N 

P 


n « N 
P P 


( 6 ) 


n 


LV. . - — in Ctj2.) 

hlo q N 

P 


n >> N 
P P 


Significant changes in high-low junction potential require that be on 
the order of N^. These simple relationships have been found to be 
satisfied very closely in the computer calculations. The detailed device 
calculations discussed in Section II show the high -low junction potential 
remaining essentially constant until high injection (n^ > lO occurs in 
the p-region. 

In summary, the computer calculations have verified and been in 
good agreement with first order models for high-low junction behavior 
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in the presence of excess minority carriers* Below high injection in 
the p-region, the high-low junction acts basically as a minority carrier 
reflecting boundary. Above high injection a part of the terminal 
voltage appears across the high-low junction. For the 100 usee lifetime, 
250 yM thick cells studied, the open circuit voltage condition occurs 
very near the high injection condition* For solar cells with longer 
lifetimes or higher p-type resistivities the importance of the high-low 
junction is greatly enhanced. Such a condition is discussed in the 


next section. 
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Figure 19. Terminal J-V characteristics for (^)' n -p with ohmic back 
contact, (b) n -p-p with 0.5 yM p region, and (c) n + -p-p 
with 5 uM p + region (no illumination). 
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Figure 20. Terminal J-V characteristics for devices shown in Figure 19 for AMO 
illumination. 
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Figure 21. 


Efficiency characteristics of the (a) n -p, (b) n -p-p with 
0.5 pM p + region, and (c) n + -p-p + with 5 yM p + region 
(AMO illumination). 
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IV. LIMITED STUDY OF THE EFFECTS OF PARAMETER VARIATIONS 
ON SOLAR CELL EFFICIENCIES 

This section, discusses some calculations performed varying various 
parameters in order to study their, influence on solar cell properties 
and efficiencies. The previous section has discussed the effects of 

'J* 

varying the back contact properties of n -p-p cells. The present 
section discusses three other types of parameter variations: (1) Lifetime 

variations, (2) Cell ' thickness variations, and (3) Cell resistivity and 
lifetime variations. As mentioned previously, all the results presented 
in this section were performed using an optical generation rate which is 
slightly (about. 10%) too large for a bare silicon surface. This should 
not effect the general comparisons presented in this section but the 
short circuit currents shown are estimated to be too large. The calculated 
curves thus lie above 10% above the expected response for bare surface 
cells and about 30% below the response for cells with optimum antireflecting 
coatings. The calculations presented in this section are for bare surface 
silicon cells with, the basic material parameters of Table I . Deviations 
from these basic parameters are discussed in each subsection. 

1. Lifetime Comparison 

The effect of p-region lifetime on the properties of the n + -p-p + 
cell was studied for values ranging from 1 ysec to 100 usee. The 
calculated terminal properties of cells are shown in Figures 22, 23, and 
24. Complete data is shown for 10 usee and 100 ysec lifetimes while 
partial data is shown in Figures 23 and 24 for a 1 useo lifetime. The 
peak efficiency is seen to range from about 7% to about 12% for lifetime 
variations from 1 ysec to 100 ysec. The increased efficiency at long 
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lifetimes results from (1) an increased short circuit current due to an 
increased carrier collection efficiency and (2) a reduced forward current 
density at long lifetimes due to decreased recombination of carriers in 
the p-type region. 

Minority carrier lifetime is one of the major device parameters 
determining the efficiency of solar cells. It is also one of the least 
accurately known of the parameters used in the solar cell calculations . 

In general it is known that .lifetime tends to decrease as the doping 
density is increased but the exact physical processes leading to such 
a decrease are not known. In order to set realistic bounds on solar 
cell efficiency, bounds on the lifetime at various doping densities 
must be established. 

The minority carrier lifetime in the n -region used in the calculations 
is 0.1 ysec. This value of lifetime corresponds to a diffusion length of 
about 10 jjM which is considerably larger than the thickness of the n + -region. 
The electric field in the n + -region due to the diffusion profile further 
increases the effective diffusion length for minority carrier collection 
at the p-n + junction. To obtain a diffusion length on the order of tfre 
n + -region thickness would require a lifetime below 1 nsec, perhaps 0.1 nsec. 
Only for such low lifetime values is lifetime in the n + -region expected to 
be an important device parameter. Computer calculations, however, need to 
be made to verify these conclusions. 

2. Cell Thickness Comparison 

Computer calculations were made varying the cell thickness to 
investigate the importance of this parameter on solar cell efficiency. 
Terminal characteristics are shown in Figures 25, 26, and 27 fop n -p-p 
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cells with thicknesses of 100 yM, 200 yM and 250 yM. Figure 25 shows 
a slight decrease in current density as the cell thickness decreases. 

This is to be expected because of the reduced recombination of carriers 
in the p-region as thickness decreases. The short circuit current is 
seen to decrease in Figure 26 as the cell thickness decreases because 
of the reduced total carrier generation rate in the thinner cells . 

The reduced forward current density with thinner cells tends' to increase 
cell efficiency while the reduced short circuit current tends to 
decrease the efficiency. The result of these two effects is a peak 
efficiency a? shown in Figure 27 which does not depend very strongly 
on cell thickness. There' is only a slight decrease in efficiency as the 
cell thickness is reduced from 250 yM to 100 pM. 

The above conclusions hold only for n + -p-p + type solar cells 
with the minority carrier reflecting back p-p + junction. For conventional 
n + -p cells the efficiency decreases much faster as the cell thickness is 
reduced below a diffusion length. The decreased efficiency in copven'tional 
cells results from an increasing normal forward current . density with" 
decreasing thickness, due to carrier flow to the back ohmic contact. 

3. Cell. Resistivity and Lifetime Variations 

This section reports on one ' computer calculation which was made 
by varying the cell resistivity and lifetime. A calculation was made 
for the 250 pM thick n -p-p cell with a p-region resistivity of 
2000 £2* cm and a lifetime of 1000 ysec. This choice of resistivity and 
lifetime was selected in order to enhance the importance of the back 
high-low junction. The values selected are approaching the highest 
resistivity and. lifetime values which can be commercially obtained with 
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bulk unprocessed silicon. It is recognized that these values may very 
well represent values which cannot presently he achieved in processed 
solar cells because of lifetime degradation during solar cell fabrication 
steps. The calculations were performed mainly t© illustrate the 
importance of the high-low junction in high-resistivity silicon. 

A comparison of the high resistivity, long lifetime device with 
a 10 ft* cm, 100 ysec device is shown in Figures 28, 29, and 30. The 
J-V characteristic of Figure 28 has several interesting features. 

Because of the high resistivity of the 2000 ft* cm cell, high injection 
occurs around 0.3 volt. This transition can be seen as a change in 
slope of the upper J-V characteristic at around 0.3 volt in Figure 28, 

The J-V characteristics for the two devices cross at about 0.45 volt 
with the 2000 ft* cm cell. having a lower current density than the 10 ftt cm 
device for large voltages. The lower curve for the 2000 ft* cm cell 
results from high injection effects and the approximate exp(qV/2kT) 
dependence at high Injection where changes in the high-low junction 
potential are occurring. The 2000 ft* cm cell also shows a slope change 
around' 0.55 volt back toward an exp(qV/kT) variation of current density. 
This arises because of current leakage through the back high-low junction 
and recombination at the back ohmic contact. Potential and electron 
carrier density plots for this device at various voltage levels are 
shown in Appendix C . 

Current density curves and efficiency curves under illumination 
are shown in Figures 29 and 30 respectively. Both devices have 
approximately the same short circuit current but the 2000 ft* cm cell 
has. a slightly larger open circuit voltage, resulting from the lower 
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current density curve of. Figure 28* The predicted efficiency of 
the 2000 ft* cm cell i$ slightly larger as indicated in Figure 30. 

The above calculations illustrate the importance of the high-low 
junction in high resistivity* long lifetime solar cells. Without 
accounting for the high injection effects at the high-low junction 
and the accompanying reduction in high-low junction potential* one 
would expect the 2000 ft* cm cell to have a lower efficiency than the 
10 ft* cm cell. This would certainly be true if the current continued 
to depend on voltage as expCqV/kT) as is observed at low voltages. 
However* the high injection effects combined with the high-low 
junction produces an exp(qV/2kT) dependence which leads to the reduced 
forward diode current -and the increased efficiency . 
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Figure 23. Terminal J-V characteristics for devices with a p -region minority carrier lifetime 
of (a) IQ ysec ( b) 100 usee and Co) 1 usee (AMO illumination). 
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Figure 24. Efficiency characteristics for devices with a p -region minority 
carrier lifetime of (a) 10 ysec, (b) 100 nsec., (c) 1 ysec. 





CURRENT DENSITY (RHP/CM 2 ) 



Figure 25. Terminal J~V character!: 
•thicknesses (a) 100 yM, 
.(no illumination). 
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Figure 27* Efficiency characteristics for (a) 100 yM ? (b) 200 yM, and 
(c) 250 yM thick solar cells (AMO illumination). 
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Figure 28. Terminal J-V characteristics for (a) 2000 1000 usee 

p-region, (b) 10 f>cm, 100 psec p-region (no illumination) 
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Figure 30. Efficiency characteristics for (a) 2000 £>cm, 1000 ysec 

p-region, (b) 10 S>cm, 100 ysec p-region (AMO illumination). 
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V. SURFACE RECOMBINATION 

Calculations have been initiated which include the effect of a 
finite surface recombination velocity at the irradiated surface. With 
this finite term the current at the surface may be expressed as follows 
(normalized form): 


d<j> 

= (p-p )S = - 2 - 
surface p 


(7) 


surface 


where S is the surface recombination velocity . Following the procedure 
illustrated in Appendix A , the variables <J> » <j> ' , and ^ may be expressed 

Jr y 

in the form 

♦pu+D * W> + s Vi)’ <8> 

where the subscript "i" refers to the iteration step. Applying this 

form to the independent variables in equation (7) one may obtain an 

expression for the correction to at the surface (5^^^ ^ ) in terms 

of the values of 4 , 4’, and at the prior iteration step. Included 

P P 

in this expression is the terra which can be expressed in terms of 

^p(lL) ^pClL 1)‘ By utilizing the "C" coefficients of Appendix A 

plus the fact that 61 ^^ = °» ^pClL-i) may be ex P ressed in tei,ms of 

6 (j)p(l L ) thus yielding the final explicit expression for The 

correction to <J> n is formed by invoking charge neutrality at the surface. 
This method provides an ongoing correction process as the device bias is 
stepped. The inclusion. of this correction in fact reduces convergence 
difficulties at the surface since the carrier concentrations at the 
surface do not need to return to equilibrium conditions. As can be seen 
from Figures H-(b) and 11(b), an infinite S demanded a large change in 


$ over an extremely small distance near the surface. 
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3 

A surface recombination velocity of 10 cm/sec was selected as a 

reasonable value and the cell described in the prior sections illustrated 

a terminal J-V characteristic (no illumination) which was within 1% of 

the prior results. However, the hole current at the n + surface was 

reduced at least two orders of magnitude. Although no calculations have 

as yet been performed under AMO illumination it can be expected that the 

2 

short circuit current may be increased by as much as 2 mA/cm . This 
estimate can be based on the consideration that the present hole current 
lost to recombination at the surface is about 2.2 mA/cm (Figure 15(b)) 
and a two order magnitude decrease of this component will add directly 
to the overall short circuit current. 
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VI . SUMMARY 

This report has summarized the work accomplished to date on the 

theoretical analysis of solar cell current -volt age characteristics 

and efficiencies. The study is only partially complete at the present 

time. The major accomplishments which have so far been achieved are: 

* - + , + + 

(1) The development of a general analysis program for n -p ana n ^p-p 
solar cells, (2) The incorporation of optical generation into the 
analysis with the capability of including antireflecting surface layers 
into the model, (3) An initial series of solar cell ' calculations 
exploring the effects of various device parameters on solar cell 
properties, and (4) The incorporation of finite surface recombination 
effects into the model. The calculations which have been performed 
illustrate the usefulness and versatility of the analysis techniques 
for studying solar cells. 

A major effort in the initial calculations has been to understand 
the influence of a back p + -p, high-low junction on solar cell properties. 
The theoretical calculations indicate that below high injection the 
high- low junction acts predominantly as a minority carrier reflecting 
boundary preventing injected carriers or optically generated carriers 
from reaching the back ohmic contact and recombining. The net effect 
on solar cell operation is a reduction in normal forward bias current, 
over n + -p cells and an increase in short circuit current. The back 
high- low junction was found to have little influence on the terminal 
properties or on efficiency when the minority carrier diffusion length 
is less than the cell thickness. When the cell resistivity is high 
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enough to. cause high injection at open circuit voltage conditions, a 
voltage decrease across the high-low junction was found to be a significant 
factor in determining the cell efficiency i To observe a significant change 
in high-low junction potential under AMO illumination was found to require 
a resistivity of .about 10 f2* cm or larger. 

Research is continuing in several areas with the present analysis 
program and future plans include the extension of the program to include 
other physical effects. The analysis is continuing into the effects of 
various parameter variations on efficiency with corrections being made 
for the slightly incorrect optical generation rate used in the early 
calculations. Calculations also need to be made exploring the differences, 
if any, when antireflecting coatings are present. The influence of 
various p-region resistivities on efficiency needs to be thoroughly 
explored since different resistivity material may lead to higher open 
circuit voltages and increased efficiencies. In this connection, 
however,. it is very important to also vary the minority carrier lifetime 
or diffusion length in a realistic manner with resistivity variations. 
Considerable thought must be given to the establishment of some upper and 
lower bounds on lifetime as a function of resistivity in order to fully 
explore this area. 

New effects which, are to be included in the analysis program are: 

(1) Heavy doping or degeneracy effects in the n + -region and (2) Two 

dimensional effects due to the geometry of the top solar cell contacts, 

The inclusion of these effects are needed to accurately describe solan 

cells. Heavy doping in the n -region is expected to be important in 

*{• 

determining the optical absorption in the n -region and the loss of- 
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carriers to recombination at the surface of the n” 1 " -region. The two 
dimensional properties are important in determining the effective 
series resistance of the cell and in determining the curve (fill) 
factor. 

In summary, the work is proceeding at present approximately on 
schedule. No major problems are evident at present which would prevent 
a successful development of a general computer analysis program for 


solar cells. 
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ABSTRACT 


A numerical technique, utilizing the computer, is developed to 
obtain solutions to the one- dimensional semiconductor device equations. 
These equations form a set of three, second order, nonlinear, coupled 
differential equations. Quasilinearization is used to convert these 
equations to a set of coupled linear equations. Iteration techniques 
are then employed to obtain numerical solutions which converge to the 
solutions of the original set of equations. These techniques h$.ve 
been applied to conventional n -n-p diodes as well as to shallow 
diffused structures of the p + -p-n + type* Some of these results are 
presented illustrating this analytical method. 



I. INTRODUCTION 


Since the origin of semiconductor device theory in the late 
1940* s, numerous attempts at solving the semiconductor device equations 
have come forth*, However, purely analytical methods have yielded 
closed-form solutions only for very simplified cases, and only over 
limited operating ranges „ Various numerical methods for the determination 
of device characteristics have been presented in the literature with the 
better ones evolving from the work of Gummel [l]* The present work 
provides accurate solutions of the device equations by means of a 
recently developed quasilinearisation process and iterative procedures . 

To utilize these techniques the one dimensional semiconductor device 
equations are normalized (as De Mari [2] and others have done) and 
reduced to a set of three coupled, non-linear partial differential 
equations . Solution of these equations for a given device structure, 
doping profile, and recombination model yields the desired device 
characteristics for either a static or transient . response . The chief 
advantage of the numerical methods presented here. over previously 
developed schemes lies in the rapid convergence of the iterations 
afforded by the quasi-linearization technique „ 

A considerable amount of work has previously been performed pn 
computer techniques for the solution of the semiconductor device 
equations (Reference E 3] presents a review of most of this work). 

Much of this work however, has been restricted to the case of zero 
recombination-generation. For the majority of devices this is not 
a. good approximation. Subsequently, one of the major objectives of 
the present work has been the development of analysis techniques which are 

7f 
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applicable to the general case of non-zero generation-recombination* This 
general case leaves one with three second order differential equations as 
opposed to the two first order and one second order* differential equations 
resulting from the zero generation-recombination model. This type of 
generalization adds another level of difficulty to the problem of 
convergence in a numerical technique [4]* Numerical solutions to the 
general case have also been reported by Fulkerson and Naussbaum [5]; 
however, their work deals with a particular device and is applicable to 
a rather restricted voltage range 0 In the present paper, the methods are 
applicable to a wide range of device structures, doping profiles, 
recombination models, external generation mechanisms, and applied voltages. 

In the analysis of semiconductor devices, the terminal currept- 
voltage characteristics are usually of major interest. The terminal 
characteristics, however, are determined by the details of the electron 
and hole densities and electric field throughout the bulk of the device 
as well as any time dependence of these quantities. Thus before the 
terminal characteristics can be obtained, the electron and hole densities 
as well as electric potential must be known throughout the device. 
Straightforward means exist for relating the terminal properties to 
the bulk carrier densities and potential. 

The basic semiconductor device equations are very nonlinear and this 
leads to difficulties- in device analysis. Even though these basic 
equations are easily formulated, they are not easily solved even with 
computer techniques. Some of the difficulties in obtaining solutions 
arise from the very rapid changes which occur in carrier densities neap 
a p-n junction. For example, across a p-n junction depletion regipn 

?0 
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the carrier densities may change by more than ten orders of magnitude 

-4 

and this may occur over a distance of only 10 cm. The total dimensions 
of a semiconductor device may, however 9 be several hundred (or more) times 
the 10~ 4 cm dimension • Thus the major changes in carrier densities and 
potential occur over distances which are usually small compared with 
the total device dimensions « When viewed on a distance scale which spans 
the dimensions of the entire device, the changes in carrier densities 
often appear as steps occurring at the interfaces between p and n regions . 
These very abrupt changes in properties at junction interfaces are very 
clearly shown in the computed curves in this paper. 

The very rapid changes in carrier densities and potential near 
interfaces provide the basis for most analytical approaches to device 
analysis. They make possible the approximate division of a semiconductor 
into space charge regions and electrically neutral (or nearly electrically 
neutral) regions. By using simplifying approximations to the device 
equations in the various regions of a device, reasonably good, first order, 
analytical approximations can be , obtained to the terminal properties of 
semiconductor devices. The abrupt changes in the character of the} device 
equations near interfaces which make possible the usual device approximations , 
give rise in exact numerical solutions to some of the major difficulties in 
solving the device equations. These problems are discussed in. more detail 
in Section III on numerical techniques. 

ft 
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II. SEMICONDUCTOR DEVICE EQUATIONS 



These equations have general three dimensional forms; however, in this 

work only the one dimensional case has been considered. In the above 

equations q is the electronic charge, e is the dielectric permittivity and 

the other terms are defined in Table 1. The net ionized impurity doping 

is represented by N(x) and can be a complicated function of x, changing 

from a positive to a negative value as an n-p interface is crossed. The 

\ first two equations (Eqs. (1) and (2)) are Poisson's equation and the 

defining relationship between electrostatic potential and. electric 

field E. Equations (3) and (4) are current density equations while Eqs. 

(5) and (6) are continuity equations for electrons and holes. Internal 

recombination or generation of electron-hole pairs due to thermal 

processes is. represented by U, while electron-hole pair generation 

due to external sources such as incident light is represented by the 

term G * 
e 
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As in most computer techniques it is convenient to work in a set of 
normalized dimensionless .variables o In the remainder of this work, 
except for explicitly indicated points , the variables will be taken 
as quantities normalized by the factors in Table 1. Utilizing this 
normalization and the interrelationship between mobility and diffusion 
constants , the device equations may be rewritten as follows : 


d 2 ip 

dx 2 


n - p + N , 


T , 3t I> 3n, 

J n * y n (n 3^ “ to? * 


(p + lE.) 
vp a- + Zx J 3 

8J 


J p = p p ^ 3x ' 3x 


_ n 

_ _ _ U _ G e - _ s 


, 3J 

iE. = _ u - g + 
at e ax 


(7) 

( 8 ) 

(9) 

( 10 ) 
( 11 ) 


Through an appropriate choice of three independent variables, Eqs. (7-11) 
can be reduced to a set of three coupled equations in three variables. 
There are several possible choices for these independent variables. 

For example ip , n, and p form one possible set. These . variables are the 
ones normally used in device analysis programs [3], The approach taken 
in this work, however, has been to introduce electron and hole quasi- 
Fermi levels, <f> n and $ , and take these, plus the electrostatic 
potential, ip, as the three independent variables. The normalized 
quasi-Fermi levels are related to the carrier concentrations as 


n = exp (iM n ) > 


( 12 ) 


and 


p - exp (<j>p-ip) . 


£ 3 - 


(13) 



6 


This approach yields three variables with approximately the same magnitude 
of variation across a device simplifying the scaling of the numerical 
techniques „ 

In terms of these three variables ip , <J> n s an <i <l> the semiconductor 
device equations (Eqs. (7)-(13)) can be written as 
2 

3 ip(x,t) _ exp ] _ exp -i|i] - N(x), (14) 

^ P 


rli 

3t 


3<p 


a 2 d> 

£ = - y (x,t)G (x,t) exp (<p -ip) + 

3x^ n e 


3<j> ^ 

n r 9ip _ n 

3x L 3 x 3x 



-Y n c*.t>uu.t), 


(15) 


and 




^ v- t> 


3 2 ^ 


3x 


■§■ = Y (x,t)G (x*t)exp(ip-4 ) + 
z p e p 


(16) 


3d> 3<p 

P r 1 


r\ 


The quantities y r and y p are the reciprocal (normalized) mobilities.. 

The nonlinearity of these equations is further compounded by the 
choice of .recombination model, character of external generation, and 
the inclusion of field and doping dependent mobilities, i_.e_. Y n and 
are functions of N(x) and 3ip/3x. For recombination, the single lev^l 
Shockly-Read-Hall model has been used. In normalized form this model 
results in the following net recombination rate: 

exp[d> -*3-1 

U _ P n 

T po [exp(\p-<t) n )-n 1 ] t T no Cexp(?) p ^)-p 1 3 



( 17 ) 
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where t , t , n„ and p., are constants depending on the location and 
po’ no 1 r l 

density of recombination levels . This has been found to be a good 
approximation to recombination in Ge and Si. For most of the numerical 
calculations , the recombination level has been assumed to exist at the 
center of the energy band gap so that n^=p^=l. This has been done for 
convenience only and the general expressions are valid for any values 
of n^ and p^. 

The mobility dependence was modeled by empirical expressions taken 
from Gwyn, et al, [6'j which are of the form 


-1 -1 2 E+C 4 1/2 

y n = U no + C 1 N / (N+C 2 ) + C 3 E E+C^ * 


(IS) 



2 E+K 4 1/2 

V [1 + V /(N+ V + K 3 E ’ 


(19) 


Normalized values of the parameters in the mobility expressions are 
given in Table II. These expressions have been found to accurately 
describe both the doping and field dependence of mobility in silicon. 
Somewhat similar expressions would also be used with other semiconductors 
with different constants. 

An external generation term G & could be developed for various 
applications; however, the mechanism studied in this work has been 
that of generation due to the full spectrum solar irradiance. This 
semi- empirical model for G^ was based on published data of absorption 
coefficients, index of refraction, and spectral irradiance. . The 
modeling of this term has been developed in work presently underway; 
however, the resulting normalized term can be expressed as follows: 


G (x) = sr / T(A)N , (A)n( X)u(A)exp[-a( A)x]dX . 

e D n. ph 

ox 


( 20 ) 



Data showing the typical dependence for silicon of G g (x) upon distance 
is shown in Fig. 1 for solar irradiation which, is present at the earth but 
outside the earth's atmosphere (so called air mass zero, AMO, condition). 
For other conditions such as absorption of light at the earth’s surface 
(AMI) or for other semiconductors, G g will be a somewhat different 
function of x» Only for monochromatic light incident on a semiconductor 
is G g (x) a simple function of distance which is easily expressed in 
functional form. In general, G e (x) is a rapidly decreasing function of 
distance from the irradiated surface , and its value is normally only 


known as a table of values. 


III. NUMERICAL TECHNIQUES 


The solution of the semiconductor device equations has been reduced 
in Section II to the solution of three coupled nonlinear equations 
(Eqs . (14)-(16)). The known values of ^ <J» n and <t> p are at the boundaries 
of the semiconductor device. Thus, the mathematical problem reduces to 
the solution of a rather complicated nonlinear two-point boundary-value 


problem. Although there are several techniques available for solving 
linear differential equations of the two point boundary value type, no 
such general techniques exist for nonlinear equations. Solutions of these 
problems generally require some type of iterative technique whereby the 
solution is approached in a series of steps. The present work utilizes 
the quasilinearization technique whereby the system of nonlinear 
equations is converted into a system of linear equations [7-9]. An 
iterative technique is then used to obtain the exact solution. The 
general principle behind this approach is expanded in Appendix A. 

As an illustrative example, the solution methods have been applied 
to the n -n-p (and p -p-n ) structure shown in Fig. 2. Since the 
electrostatic potential (i|0 is arbitrary to within a constant, it may 
be taken as zero at x = 0. Its value at x = L is then the built-in 
voltage plus the applied terminal voltage. The boundary conditions of 
the quasi-Fermi levels and $ are then determined from the defining 
Eqs. '12 ) and (.13) and the values of the carrier concentrations, n and p, 
at the boundaries. If ohmic contacts are assumed at both end points, 
the carrier concentrations are fixed at their equilibrium values which 
are independent of the applied voltage. This leads to the condition 

that <b and <b are both constant at x = 0 and increase at x = L 

T n p 
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directly with the applied terminal voltage. The restriction of ohmic? 
contacts at x = 0 and x - L may be loosened to allow for a finite 
surface recombination velocity at this end point. However, to more 
clearly illustrate the analysis technique, ohmic contacts will be 
assumed. 

To proceed with the technique, the coupled device equations can be 
written in the following functional form: 


3 ip(x,t ) _ f . , , x 

— - — £ — = F.(ift,<fc ,6 ) 


3 6 (x,t) 


= f 2 U» * n . 



F3U9 * n » $ p » 


(23) 


The functions F , F 0 , and F correspond to the right hand sides of 
Eqs. (14), (15), and (16) respectively. 

In the quasilinearization technique (as discussed in Appendix A) 
the nonlinear differential equations , instead of being solved directly , 
are solved recursively by a series of linear differential equations. 

The linear equations to be solved are obtained from the first order 
terms of a Taylor's series expansion of the nonlinear part of the 
original differential equation (terms F^, , and F^ above). This 

technique is a generalized Newton-Raphson formula. for functional 
equations [7-9]. In order to illustrate the quasilinearization technique 
as applied to the semiconductor device equations, Eq. (22) is considered 
as a typical equation. 
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Applying the quasi linearization technique to Eq. (22) gives 

_ 3(j> . . 3% 

|- m+1 _ n ,m+I -. T n ,m+l _ . 

L 8t 3t JY n " 9x 2 2|m 


2 2 1 

+ 3ip ^nr^m+l + &4> ^n,m+l ^n,m^ 


SF 9F 

+ 3A ^m^p.iiw-1 ^p + 1 ^ m ^m+1 


+ ^L ~ K,J- (24) 

» 

where the subscripts m and m+1 indicate the m and m+1 iteration steps . 

Approximate solutions to to , to , and to are assumed known in order 
" r m Y n,m p >m 

to generate the solution at the m+1 iteration step . The partial derivative 

terms of F 2 such as 9F 2 /3to n l m are also assumed known since they are 

evaluated using to , to » and to . Also to* and to' represent first 
6 r irT y n,nf p s m r Y n K 

derivatives with respect to distance.. Equations similar to Eq. (24) can 

easily be obtained for the other two equations (Eqs. (21) and (23)), 

By the above technique the three device equations are converted into 

three linear differential equations. Starting with some initial approximation 

to to, to , and to at iteration step m = 1 the functions to, , to ,, and to , are 

n’ T p r r l’ n,l’ pjl 

known , and to 2 » 2 J an< ^ ^p 2 are calc ulated. These new solutions are used to 
generate a better solution set with the process continuing until the 
solutions converge within some specified degree of accuracy. In solving 

% 

the linear differential equations in an iterative manner, it has been found 
convenient to express the equations in terms of correction terms at 


the m'th iteration step as 
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6 *m = *ml ' V 


66 — 6 _ — 4* ^ 
Y n 5 m+1 n 3 m 


§6 =6 — 6 fl 
p,m T p,m+1 p,m 

Since the solution is numerical* discrete points in space and time must 
be used at which the values of ip 9 and 4>p are calculated * These 
points are identified as i*k; where i is an integer specifying the 
x-coordinate and k is an integer specifying the time coordinate* Applying 
these definitions to Eq e (24), it may be rewritten as 


36 d> 3 2 6cb 

r®£i JL l y - _E = h 

L 9t . , 3t . , JT n . _ 2 . , 1 

m,i,k m,i # k m,i,k dx m,i,k 


2 • v 

* m,i *k 


'm,i,k 


a *\ ) + f2 

Y " 3t m,i ,k 8x 2 m,i,k ' ** 


m,i ,k 


6i(j . , 
m,i ,k 


— 64> + —• H 

M 15 ♦ 1 3^ p ■> L, 

n n ,i,k ”’ 1,k P m,i,k "’ i ’ k 


df + »f r & K 

m,x,k m,x 5 k T n'm*i,k m*x,k 


Similar differential equations may be obtained for dip and.6$ . The 
solution to these equations give the corrections to ip , 4> n and $ at each 
iteration step. 

The system of differential equations has been solved by the finite 
difference approach using approximations of the Crank-Nicolson type [101 
for the time dependent . case and of the three point difference type for 
the time independent case. Following the use of finite differences to 



approximate the derivatives, a set of three coupled finite difference 
equations results, which are of the form, 
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6 *i,i - s h + 4 v + 4 5 v + 4. + 


>n. , = A 21 S *i + A 22 S V + A 23 S *p- » A 24 + 
X-l 1 *X 


B 21 6 *i+l + B 22 6 *n. n > 

x+1 


6 *p. n = A 31 6 *i + A 32 6 *h. + A 33 6 *p. + A 34 + 

*x-l X -t x 


B 31 + s 33 s +p. • 

r x+l 

The m and k subscripts indicating the iteration step and the time step 
have been omitted for clarity since it Is understood that the correction 
terms are being evaluated for the same time and . iteration step. The 
and B^_ coefficients are complex functions of position and of the 
potentials (>|>, <t> n> and at previous time steps and previous iteration 
steps. The above equations relate the corrections to tp, and at 
spatial point i to the corrections at spatial points i-1 and i+1. The 
A and B coefficients are slightly different for the time dependent and 
time independent cases. Exact expressions for the time independent case 
are given in Appendix B. Expressions have also been derived for the 
time dependent . case [11]. 

The application of the above equations to a p-n junction is straight- 
forward in principle. However, there are several difficulties which must 
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be overcome before solutions which converge can be obtained . One of the 

most important of these problems concerns the selection of spatial points 

at which ip, <p , and d> are calculated. In a p-n junction at low forward 
Y n p 

bias, most of the changes in ip, 3bd <p^ occur over very small distances, 

either within the depletion region or for small distances on either side 
of the depletion region. To achieve accuracy in the numerical calculations, 
any changes in ip, <p , or from one spatial point to another spatial point 
must be small. This requires that very small steps be employed within 
space charge regions. On the other hand outside these regions the quantities 
change much less rapidly and efficient utilization of computer time requires 
that much larger step sizes be employed. It is necessary then to use a very 
nonuniform spatial step distribution . Because of this difficulty Fulkerson 
and Nussbaum [5] in their numerical calculations on a p-n junction divided 
the device into three regions and used the Debye length as a spatial 
normalization constant within the depletion region and the diffusion length 
as a spatial normalization constant outside the depletion region. The 
approach used in the present work has been to use the same normalized 
equations throughout, but to generate a nonuniform spatial step distribution. 

The generated spatial step size distribution is such that small 
step sizes are taken within regions of space where ip, ip^, or <p^ are 
changing rapidly and larger step sizes are taken when these functions 
change slowly with position. This step size selection process is similar 
to that developed by De Mari -[2]. Since some of the final results (i_.je. 

and . J ) depend upon an integration of the exponential functions of 
ip, <p n » and ipp, the criteria on step size is selected so as to maintain a 
constant relative error within similar numerical integrations. The 
relative error is defined as the fractional difference between a 
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trapezoidal and a parabolic numerical integration of exponential functions 
of these variables, and can be set to some convenient value. This in turn 

can be related to the step size through the radius of curvature of. the 

♦ 

exponential function [2,11]. At the initialization of an analysis, first 

order approximations are used for ip, and $ within space charge regions 

while final data at a prior voltage increment is used' in subsequent voltage 

steps. This allows the step selection to adjust automatically to the 

changing situation as the applied voltage is incremented. The step 

selection has been found to be one of the crucial factors in obtaining 

proper convergence. This becomes particularly apparent in very shallow 

diffused devices such as solar cells, since the potentials vary through 

a considerable range within very narrow regions. On the other hand there 

are basic computer limitations due to round off errors on the minimum 

size of a step. Thus a, proper balance of step size and number of steps 

in a particular region must be maintained. For a typical device run, 

calculations can be normally made with 1000115 steps across a device, 

with the normalized step sizes ranging from around 10 ^ near the center 

-15 

of the device to around 10 near a surface of a shallow diffusion. 

Double, precision is also used in the calculations to extend the range 
of allowed step sizes. 

The coupled difference equations for ip, if> n and ere tridiagonal 
matrices. A special technique has been developed for solving a 
tridiagonal matrix [9,12] and in the present work this technique has 
been extended to the solution of coupled tridiagonal matrices [11]. 

For given boundary conditions on ip, and $ the solution of the 
coupled difference equation is a straightforward procedure. In the 
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present work the boundary conditions are fiifi = 5$^ = 6$ = 0 at both 

boundaries of the device since the potentials do not change from one 
iteration step to the next* 

With the aid of the flow diagram in Fig. 3, the general procedure 

for the analysis, of a p-n junction such as shown in Fig. 2 can be 

summarized. First, information on the device structure such as region 

widths, doping levels, applied voltage, etc., is used to calculate 

first order approximations to the potential and quasi-Fermi levels 

throughout the device. For zero applied diode voltage, the quasi-Fermi 

levels are constant across the device and the first order approximation 

to if> is obtained from standard abrupt junction or linear graded junction 

theory depending on whether the junction is abrupt or diffused [13]. 

For applied voltages larger than zero , the first order approximations 

are taken as the converged solutions at the previously calculated voltage 

point with the voltage increment linearly added to the previously 

calculated potentials. The initial approximations are used to develop 

a step distribution and to obtain values for the A and B coefficients 

formulated in Appendix B. With these coefficients the linear difference 

equations are used to calculate corrections to iJj, <f» n> and <|> . Thege 

corrections are applied to the old values of Uj , <J> , and $ and the net 

Ti P 

result continues in the iteration loop to form the basis for a new 

calculation of the A and B coefficients from which more accurate values 

of iff, 4> n , and are obtained. The process is repeated until the 

largest correction on either if;, <f> n , or $ is less than some previously 

-4 

set quantity (usually 10 in this work). When the largest correction 


at any point approaches zero the solutions satisfy the original semi- 
conductor device equations . With ib , <j> , and <j> known , the carrier and 

n p 

current densities can be evaluated. 
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The solution for any number of voltage values can be obtained in a 
single computer run. The number of iterations required at each voltage 
step is around five when the voltage is stepped in increments of 0.05 
volts. A complete forward I-V characteristic can be obtained from 
0.0 volts to 0.7 volts in .05 volt increments in about three minutes 
of total computer time on an IBM 360 computer. Successful runs have 
been made up to 2 volts applied forward bias; however, the major 
interest of the present work has been at voltages up to 0.7 volts. 

A flow diagram for the time dependent calculations is shown in 
Fig. 4. The calculations for this case are carried out in a similar 
manner except that each calculation is for a different point in time. 

The time increments in the program are similar to the voltage increments 
in the time independent analysis program. Most of the work has been 
with the time independent program; however, sufficient calculations 
have been made to verify that the techniques can be successfully applied 
to time dependent problems . 
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IV. EXAMPLES OF COMPUTED RESULTS 

The methods outlined in the previous sections have been applied to 
a wide range of .abrupt junction n + -n-p* structures with various doping 
densities and region widths* Some of the calculated results have been 
previously reported [14]* As illustrations of the results obtainable 
from the program, one series of those reported results will be reviewed 
and more recent work involving p -p-n shallow diffused devices will 
be presented. These were selected in order to illustrate the versatility 
of these numerical methods over a rather wide range of diode structures. 

Calculated data for two types of devices is included in this work. 
These are, first, an n -n-p diode with a lightly doped n region and 
wide n + and p + regions. In this device the impurity doping (N(x)) is 
assumed to abruptly change values at the junction interfaces . Other 
pertinent data on the device is shown in Table III . For this n + -n-p + 
device the center n region is about a minority carrier diffusion length 
in width while the n + and p + regions are considerably wider than a 
diffusion length. Calculated plots of potential and quasi-Fermi levqls 
are shown in Figs. 5 and 6. The rapid changes in potential at the 
junction interfaces can be seen in Fig. 5, where the potential change 
appears as a step in the curves. Due to the scale of presentation the 
space charge regions appear as steps in the figures even though each 
region has a substantial number of computational points which can be 
retrieved via a numerical printout. Although the high injection region 
was reached as evidenced by the rise in majority concentration of the 
n type region seen in Fig. 7 S the high bias condition was mostly 
dominated by series resistance as seen from the linear 

\! 

i 
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increase of potential within the center region at 0,8 volt bias (Fig. 5). 

The carrier density curves shown in Fig. 7 are essentially as one wou^d 
expect from diffusion model calculations ? Details of the changes ij} 
carrier densities within the space charge regions can again not be 
seen due to the scale of representation 9 so this region corresponds to 
the essentially vertical steps in the carrier densities seen in the 
figure. The electron and hole quasi-Fermi levels as shown in Fig. 6 are 
essentially as one would expect from first order device models. As shown 
in the figure, the quasi-Fermi levels are essentially constant within the 
depletion region. 

The calculated results for a shallow diffused p -p-n device are 
illustrated in Figs . 8-12 . The analysis of this device has been under- 
taken as part of a study of silicon solar cells which is presently in 
progress [15]. The dimensions and doping densities correspond closely 
to typical solar cell structures . This device has been studied both 
without and with optical generation of electron-hole pairs. The potential 
and quasi-Fermi level plots at various voltage levels for this device are 
shown in Figs. 8 and 10. Again the changes at the junction interfaces 
are very abrupt on these scales and appear essentially as vertical lines. 

The n + region is so narrow that no detail of this region can be seen in 
either Fig. 8 or 10. An expanded view of the potential variations with 
applied voltage around the p-n + junction is shown in Fig, 9. Below high 
injection in the center region, which occurs at around 0.6 volt, essentially 
all of the potential change with increased voltage occurs across the 
p-n + junction. For voltages above 0,6 volts approximately half of the 
increments in terminal voltage appear across the p + -p junction, as 
predicted by first order models of diode behavior. 




Electron and hole densities are shown in Fig. 11 for the entire 
device, and an expanded view at the p-n junction is shown in Fig. 12. 
The boundary condition of an ohmic contact at the surface forces the 
surface values of p to be independent of applied bias. Expanded views 
are quite necessary in this device to see any variation of n and p at 

t 

the p-n junction interface. 

This same device was analyzed with the inclusion of a generation 
rate due to solar irradiance. Fig. 13 illustrates the overall change 
in the carrier distributions and Fig. 14 expands this data about the 
n + -p junction near the irradiated surface. These figures can be 
compared directly with Figs. 11 and 12. The inclusion of a generation 
term represents a rather heavy perturbation of the program but it was 
found quite easy to add the generation rate in steps to an already 
convergent solution at some high forward bias . After stepping in 
powers of 2 from Ge/128 to Ge/1 at the same forward bias, the voltage 
could then be incremented downward in the same fashion as before in 
order to obtain a full characteristic. Once the program converged with 
the full generation term, it incremented in voltage to the short circuit 
condition without any difficulty. 

Although no results are shown in this paper, the programs also 
evaluate the terminal current density at a given voltage, and generates 
terminal I-V characteristics. The evaluation of current density is 
straightforward from the device equations after ^ and <fi are known . 



V. SUMMARY 


The technique of quasilinearization and the numerical techniques 
discussed herein have been found to be a very good approach to 
semiconductor device analysis . The techniques have been found to lead 
to rapidly converging solutions under a wide variety of physical 
conditions . These cover the range of devices which are narrow compared 
to the minority carrier diffusion length to devices which are very 
wide compared to the minority carrier diffusion length. The approach is 
also able to very easily handle high injection effects in p-n junctions, 
high-low junctions, and optical generation. The exact analysis is 
providing very useful information on the operation of solar cells. 

With the analysis program a much more complete understanding is possible 
of solar cell operation and especially of the interactions between p-n 
junctions and high-low junctions in solar cell . structures [15], 

The results reported herein are but a small sample of the computed 
results which have been obtained on semiconductor devices. The included 
results illustrate the application of the program to two physically 
very different diode structures. The analysis program is able to analyze 
both types of devices with only minor changes in input parameters. 
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Table 1 Normalization table 


Normalized 

Quantity 

Description 

Normalization 

Factor 

X 

position coordinate 

L d < = /eV t /qn.) 

t 

time coordinate 

L >o 

T no * T po 

lifetimes 

L >0 

♦ 

electrostatic potential 

V (= kT/q) 

6 3 0 
r n p 

quasi-Fermi levels 

h 

V A 

applied (or terminal) voltage 

V t 

E 

electric field 

V t /L D . 

n »P 

carrier densities 

n„ 

i 

N, N d , N a 

net impurity 5 donor * and acceptor 
densities 

n. 

i 

J ’ J n* J p 

total, electron, and hole current 
densities 

-’VA 

D , D 
n* P 

carrier diffusion constants 

2 

D (= 1 cm /sec 
o 

JJ , u 
n p 

carrier mobilities 

V v t 

U 

generation -recombination rate 

Vi /L D 

G 

e 

external generation rate 

D oV L D 


/#/ 
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Table II . Constants 

for mobility 

expressions 

(normalized values) 

C 1 

81 

K i 

350 

C 2 

2c 16x10 8 

K 2 

7.0xl0 8 


— 6 


_ -6 

C 3 

1.76x10 

K 3 

1.93x10 

V 

l«70xl0 4 

K 4 

1.76xl0 4 

C 5 

1.25xl0 4 

K 5 

1.18xl0 4 

-1 

y 

34.93 


12.42 

no 


P° 
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Table III. Data on two diode types studied 


+ + 
n -n-p 

Device 


Device 

+ + 
p -p-n 

(Solar Cell) 

n + doping 

10 18 /cm 3 

p + doping 

10 18 /cm 3 

n doping 

2 . 5xl0 15 /cm 3 

p doping 

1.3xl0 15 /cm 3 

p + doping 

10 18 /cm 3 

n + doping* 

10 20 /cm 3 

n width 

250 ym 

p + width 

0,5 p 

n width 

100 ym 

p width 

249,25 ym 

p + width 

250 ym 

n + width 

0.25 ym 

T 

po 

, -7 

10 sec 

T 

po 

-7 

10 sec 

T 

no 

—6 

10 sec 

T 

no 

-4 

10 sec 

*n region is a Gaussian diffused 
the density at the surface. 

region . The value 

of doping density is 
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APPENDIX A 

The quasilinearization or Newton-Raphson approach to nonlinear 
problems was first introduced by Bellman [9]) extended by Kalaba [10] 
and applied to many engineering problems by Lee [113 » To Illustrate 
the mechanics of the technique, consider the equation 

L[X(t)3 = g(X,t) (A-l) 

where L is a linear operator (partial or ordinary) operating on the 


desired function X(t) and g is some nonlinear function of X. Performing 
a Taylor series expansion of the nonlinear function about a given 
initial estimate X^Ct), one obtains 

L[X(t)3 = g(X o ,t) + (X-X q ) g'(X o ,t) + CA-2) 

where the prime indicates differentiation with respect to X. Neglecting 
second and higher order terms of the expansion, Eq, (A-2) becomes 

L[X(t)3 = g(X o ,t) t (X-X o ) g‘(X 0 ,t), (A-3) 

which is a linear equation which can be solved for X(t) = X^(t). Using 
this result to formulate 

L[X 2 (t)3 = g(X 1 ,t) + (X 2 -X 1 ) g , (X 1 ,t), CA-4) 


the iterative process can be continued until the desired solution X(t) 
is obtained. Of course, this procedure assumes convergence of the 
sequence of functions fx }. Although at first this seems to be a rather 
severe restriction , for many of the commonly occurring linear operators , 
the sequence does indeed converge. In fact, for all functions attempted 
by the authors, convergence was accomplished. It should also be noted 
that if convergence does occur, then this convergence is quadratic in 
nature [9,11], That is, as convergence is approached, the maximum error 



(or correction) in a given iteration is proportional to the square of 
the maximum error of the previous iteration. 

In summary, the quasi linearization technique is essentially a 
Taylor Series expansion (in function space) of the nonlinear portion 
about each of the individual functions making up the nonlinearity . 

The resulting linear equation is then solved in an iterative fashion 
subject to the boundary conditions imposed by physical considerations. 
As evidenced by the equations in the body of this paper, the 
quasi linearization technique is applicable to multi-dimensional and 
multi-functional nonlinear equations. 
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APPENDIX B 


Expressions for the A and B Coefficients for the Time Independent Case. 


a. , s. . 3F. 

A?-, = [-—■ + ( s .+s . ,) 


11 a . 

i 


2 i i-1 3^ 


], 


(B-l) 


s. . 3F 

A 1 = ~ — (s +S ) — 

A 12 2 K i i-l ; 3* 

n 


(B-2) 


s. . 3F. 

A 1 - (s +s ) — — 

A 13 2 ^ i + i-l J 34 


(B-3) 


A 14 = 


(B-4> 


B 11 * ‘^V 


(B-5) 


A 21 = ( 


s. 

X 


3F, 


2+s . 


i 3<t>' 
Y n 


s. + s. . 3F o1 

.) J X Xz l _ g. 

; 1 a. 3ip I 
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a. . 

i-I 


i 3F 2 

A*,) 


a 


2 ir 3^» 
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A 22 - 


s. 
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o Hi 

2+S i 3t|> ' 
Y n 


-) { 


s.+s, n 3F 0 
x x-1 2 


a. 
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3<J> 
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2 ( a « + 1 ) 

+ — ■ — + 

s.ou 

X X 


2 Scf. T 

a. 

x n 


- A 


i 3F 2 


12 Sip* 
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s i (s i +s i-l ) r 3F 2 


3F„ 


2 + s. 


i 30' 
n 


s. SF 
x 1 

2 3 4 1 


3F„ 


30' 


} . 


C B— 8 ) 


s.(s.+s. , ) t, s. 3F„ 

1 1 1-1 r~_ ' 1 - 1 i-_ ...1 l 


24 3F 2 

2 + s i ai 7- 
n 


< CF 2-*n ,;l i -4 tr i'*" 3 


30* 
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3F, 
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i 30’ 
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i 3$ ' 
P 


s.+s. 3F- 

.)(, 1 . 1 ~ 1 —1 

>X a. 3iJ> 


.—■■■■. ■) £ 

V 2 J 
a, 
x 


- A 


i 3F 3 


11 3t p' 


i , 


CB-12) 


32 


s. _(s.+s. ..) 

( 1 - 1 \A X ■ ) { 


SF, 


2 + s. 


i 34)’ 
P 


3F, 

v 

30~ 


s « 3F 

3F 0 

x 1 

3 

2 34 

Y n 

30* 

X 

i 


] , (B-13) 


33 * ( ‘ 


9F 3 
2 + s i 8V 


(Vi,!!! 
1 2 ’ >*• 
“i P 


s. + s. - 3F- 

> i- — ^T 

1 a . 30 

x r p 


a >s . 
x x 


- A 


i ^3 
13 30 


} , 


(B-14) 




30 


*34 - ( 


S, _(s Q +So 
1-1 1 1- 

-1 

3F 3 


2 + S 1 34 1 * 
P 

i 


-) 1 [r 3-+'p J . - r CF r*" ] . 

* T 1 


' F 3i 

34 * 


. a.. , 

t ,1 _ r t ii-l\ 3 

B 31 ‘ Cs i-1 ( a. } 34' 


3F 3 

1 / 12 + & iW 

P 


} . 

(B-15) 

(B-16) 


. s. . 

n 1 - r 1-1 
33 “ L 3F„ 


3 [■ 
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In the above equations s. is the magnitude of the step In going from 
x. to x. , and a. is the ratio of s. to s. , . Expressions for the 

i i+l i i i-l 

various partial derivatives of F^ , , and F^ are given in the following 

tables : 

Table B.l Partial derivatives of F^ 


F^*, 4> n » 

4p) = expC4)-4> n 3 - exp[ 4>p~40 - N(x) 

3F 


1 

34 

= expCip'^J + exp[4p~40 

3F 


1 

3 <fa 
r n 

= -expC4)-4> n 3 


3 F 

— - = -exp [ 4-40 

34 p P 


/<?2 
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Table B . 2 Partial derivatives of F^ 


d<J> dy j . d<p 

'IsrViT'^df* Y n (K)G e (x) e xp(t n -*) 


exp[<p — ] - exp[0 ~\p] 

+ y (x) ^ 

‘ l t [l+exp(ip-<p )]+t [l+exp(<j> -tp)3 


Po 


n no 


(t +t )[exp(<f> -ip) - exp(<p -(p)] + 2 x [exp (4 -4 )-l] 
o no po P n po P n 

- y (x) — — 5 

^ n {t [l+exp(ip-<p )] + t [l+exp(<p -ip)^} 

po n no P 


-y (x)G (x)exp(<j> -ip) 
n 0 n 


(t + t )exp(<p -ip) + t exp((j> +<j> -2*p) + x E2-exp(<j> -<p )] 
3F no po r r n r no r T n T p po p n 

w - \ U) 

T n 


{x^Cl+expC ip-ip^ ) ] + X nQ El+exp( <P n -ip)]} 


po 


+Y n ( x) G e ( x) exp ( 4> n -ip ) 


g F 1 (T no +T po )eXp( V^ + V^V^n 5 + T no exp( ^ n + ^~ 2 ^ 

-r-7- = - y (x) — -z 

3 *p fT po [l+exp(<p-.j> n )] + x no [l+exp((p p -ip)]} 


9F 2 d^Cx) 
3tp' _ dx 


3F 0 n dy (x) d<p (x) 

2 _ X n ry n 

3 d> 1 - y ( x) dx dx 

n n 


dipCx) 

dx 
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Table B.3 Partial derivatives of F, 


dy dij> d<fi 

F„(tl),ib ,<b ,Ui’ ,4 ' ) = + — j-E* " “Y ( X ) G (x)exp(t|i-$ ) 

3 y,v n* r p r y p dx dx dx dx 'p e c p 


exp(ip-<(> n ) - exp(^-ifip) 

+ V X) i po l- 1+e *P^-* n )] + t no Ll+expC(() p -^). 


(t +t )[exp(i/>-4 ) - exp(i*)-4 )] + 2x [exp(4 -6 )-l] 
3F no po n p no P n 

-TV- = Y_,(x) — 2 

^ P {t [l+exp(t*>-if> )] + x [l+exp(<}> -i*)]}" 

1 po n no P 


-Y ( x) G e < x) exp ( ) 

(t +t )exp( 1*1-4 ) t x exp{4 -4) + t exp(2i|i-4 -<#» ) 
dfo no po n no P n po n P 

3 = -y ( x) 5“ 

d *n P {x [l+exp(i*i-4 )] + x [l+exp(ij> -ip)]} 

po n no P 


(x +x )exp(i*>-<|> ) + t expC 2i|t~d ) + x [2— escp< 4> — <f> )] 


no po ‘ P po ' ' P no 

tt— - y (x) 7T 
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n 'p no 
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FIGURE CAPTIONS 


Spatial dependence of electron-hole pair generation rate for 
full spectrum (AMO) irradiance on silicon «> 

The basic n + -n-p + (or p + -p-n + ) diode structure. 

Flow diagram of semiconductor device analysis program 
(Time independent case). 

Flow diagram of semiconductor device analysis program 
(Time dependent case). 

Electrostatic potential variation in n -n-p device. 

Quasi-Fermi potential variation in n -n-p device 
(measured above equilibrium value of -Q.466v). 

Hole and electron density variation in n + -n-p + device. 

Electrostatic potential variation in p + -p-n + device. 

Electrostatic potential variation about n -p junction. 

< 4" -j- 

Quasi-Fermi potential variation in p -p-n device (measured 
below equilibrium value of 0.466v). 

+ + 

Hole and electron density variation in p -p-n device. 

Hole and electron density variation about n + -p junction. 

Hole and electron density variation in p + -p-n + device 

under full spectral irradiance (Irradiated surface at 
x = 250 microns). 

Hole and electron concentration about n + -p junction near 
irradiated surface. 
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Figure 2. The basic n -n-p Cor p -p-n ) diode structure. 
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Figure 3 „ Flow diagram of semiconductor device analysis program 
(Time independent case). 








Figure 4. Flow diagram of semiconductor device analysis program 
(Time dependent case). 
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Figure 6. Quasi-Fermi potential variation in n -n-p device (measured 
above equilibrium value of -0.466v). 
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Appendix B: Optical Carrier Generation in Silicon 

P. M. Dunbar and J. R. Hauser 
Electrical Engineering Department 
North Carolina State University, Raleigh, N. C. 27607 

ABSTRACT 

A numerical analysis ©f the spatial dependence of optical 
generation of electron-hole pairs in silicon due to full solar 
spectral irradiance is reported. The analysis includes the 
utilization of published experimental data on solar spectral 
irradiance, index of refraction, and absorption coefficients. Of 
particular interest is the illustration of the depth to which one 
must collect optical carriers in order to utilize given amounts of 
the incident photons. A determination of reflection and transmission, 
coefficients is made for various ' thicknesses of SiO and Si0 2 
antireflection coatings and optimum thicknesses for photon 
transmission are reported. Absorption in the SiO layer is included. 



I . GENERAL FINDINGS 


Electron hole pairs generated in a semiconductor arise from the 
interaction of photons with the valence band electrons in the semicon- 
ductor, The irradiated material is .characterized by the absorption 
coefficient, o, which is a rather complex function of the incident, 
photon energy. For full solar spectral irradiance, the number of 
incident photons as well as the reflectance properties of silicon are 
quite complex functions of the photon energy. Also at the higher 
energies of the solar spectrum the internal quantum efficiency exceeds 
unity. For these reasons an analytical development of the optical 
generation rate due to full spectral irradiance has not to date been 
undertaken. The present work incorporates experimental or empirical 
data concerning each of the above parameters in order to formulate the 
spatial dependence of the optical generation rate.. This is performed 
for both air mass one (AMI) and air mass zero (AMO) solar spectral 
irradiance. Also calculated are the effects of both SiO and Si0 2 
antireflection coatings on generation rate and the optimum thickness of 
these layers for maximum photon transmission into the silicon substrate. 
Absorption in the SiO film is included. Prior work in this area, has 
been limited to equal energy spectral response and does not ■ integrate 
into the response the full spectrum of solar irradiance, reflection 
losses, quantum efficiencies, or absorption in the antireflection 
film [1-8]. 

For monochromatic radiation an incident wave will decay in energy 
exponentially with respect to distance . This 'can be expressed as : 
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P(x) = P(0) e~ aX , (1) 

4irk 

where a = — ~ is the absorption coefficient , and k , the extinction 
A c 

coefficient, is the imaginary part of the index of refraction. 

Neglecting internal reflections the maximum current resulting from 
this incident energy may be expressed as : 

J Qe = / / qT( X ) n A £l-exp(-aU)x)] dA dx, (2) 

where T is the transmission coefficient, is the incident 

photon density per wavelength interval , and is the internal quantum 

efficiency. The generation rate is also expressed as 

G (x,A) = T( A) N fX) n, a( A ) exp(-a(A)x). (3) 

e phi A 

A detailed discussion of the methods used to model the various components 

in the above expression is given in the next section. Figure 1 however, 

illustrates the resulting spatial variation of the calculated generation 

rate for a bare Silicon surface. Shown in Fig. 2 is the effect of 

SiO and SiO^ single layer anti-reflection coatings under AMO conditions. 

22 3 

For AMO conditions the generation rate is seen to exceed 10 /cm *sec at 
the surface and to decrease rapidly with distance into the semiconductor. 

A significant decrease is seen in generation rate as close as 10 cm (10 pm) 
from the surface. 

A graphic illustration of the very rapid decay of the generation term 
can be obtained through an integration of the generation term versus 
distance. This calculation yields the maximum optical generation 
current as a function of the device depth over which optical carriers 
are collected. The result is presented in Fig. 3 in terms of the 
collection efficiency as a function of device thickness. . 
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It is assumed in this calculation, as in the others, that the back contact 
is transparent. The total optical current In this work is defined as the 
total photon current which Is present in the visible spectrum above the 
silicon band edge Cl. 11 eV in this work). The results for the other 
geometries and AMI irradiance are within 2 or 3 per cent of this curve . 

Table 1 summarizes the various . optical parameters used in the 
calculations and results obtained. The total surface loss for a. particular 
condition is the percentage of incident photons which are not transmitted 
into the device. This loss includes reflection and any absorption in the 
antireflection layer. Since the reflecting nature of silicon and eny 
antireflection layer is strongly dependent of incident wavelengths the 
total surface loss will vary from AMI to AMO due to the difference in the 
two spectral densities. Results similar to these were presented by 
Milnes and Feucht E93; however, their work assumed the SiQ layer to be 
fully transparent, while in this work losses in the SiO layer have been 
taken Into account. The theoretical calculations predict an optimum 

Q O 

antireflection coating of about 80QA for SiO and about 110QA for SiOg. 

The predicted surface loss for SiO layers is slightly less than that 
for Si0 2 . 




II. REVIEW OF FACTORS DETERMINING GENERATION RATE 


It is known that monochromatic photon absorption has an exponential 
dependence upon the absorption coefficient* a. Although it is possible 
to. determine a functional dependence of the absorption coefficient upon 
photon energy it is quite complex, so experimental data on a vs . X 
has been used. This data is that of H. R. Phillip [10], and Dash and 
Neuman [11]. However, in using this data it is necessary to review 
the major absorption processes in order to justify any assumptions 
pertaining to quantum efficiency. The single most important process 
for solar cell operation is that of intrinsic absorption, which is the 
generation of hole-electron pairs by photons of energy larger than the 
band gap energy . This is the most predominant mechanism in the solar 
spectrum range of wavelengths . 

At energies above about 3/2 the band gap energy there is the 
possibility that hot carriers may create secondary pairs via impact 
ionization and thus raise the quantum yield above unity. This becomes 
appreciable in silicon for incident energies above 3 ev and is included 
in this present work by the following empirical model of Shockley [123: 

= 3-2 exp(~ "-| — g^0 , hv > 3.3 ev (4) 

Although this increase in quantum yield is included, it can be considered 
somewhat secondary due to the low spectral power above 3.3 ev . Also the 
high absorption coefficients in this range of energy imply high surface 
absorption and probable subsequent loss due to surface recombination which 
would b© present « 
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At energies at and below the absorption edge there are several 
mechanisms. The more important of these is free carrier absorption. 

This absorption is seen as a rise in a for photon energies near and 
below band edge. Of course, this mechanism has a quantum yield of 
zero. This zero yield is simulated by omitting any appreciable 
absorption coefficient data for energies below the band edge. However, 
free carrier absorption does still occur for energies above the band 
edge and could conceivably lower the effective quantum efficiency . 

Data by Spitzer and Fan [13] indicates a dependence of a on the square 
of the wavelength and on the carrier concentration. An extrapolation 
of their data to energies above the silicon band edge indicates that 
a is at most 50 cm” 1 . Thus, the typical device will be effectively 
transparent to this effect. 

Other low energy absorption mechanisms include lattice vibrations, 
exciton generation, band tailing, isoelectronic traps, and extrinsic 
centers (non ionized impurity centers) [14,15]. However, all of these 
are appreciable only into the infrared and even then have extremely 
low absorption coefficients. These were neglected in the development, 
of the generation term in the same manner described above. 

Reflection from the surface of the sample must be taken into 
account, A polished silicon surface, calculated as shown below, 
will lose in the neighborhood of 35% of the incident photons due to 
reflection. With a single layer antireflection coating this loss is 
reduced appreciably as shown in Table 1. For the structure illustrated 
in Fig. 4 the optical coefficients can be written as: 
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R = 


r i2 2 + 2r 1 ^ 23 exp(-2e)Cos[26- (^ 23 - ex P ( ~ 4&) 

1 + ^ r 12 r 23 )2eXp(: ” 4e) + 2r i2 r 23 eiip(_26)COSt - 25 "^23 + *12^ 


(5) 


T = 


( t 12 t 23 ) exp(-26) 


1 + ( r ’ l2 r 23 )2exp ^ -4e) + 2r ' 12 r 23 eXp( ‘" 2 ^ CoSt2S ” ( ’ <J '23 + *12^ 


( 6 ) 


T + R + A = 1, 


(7) 


12 


^l 2 - ^ - k c/ + (2k c2 n l^ 
(n 1+ n 2 ) 2 + k c2 2 


1/2 


( 6 ) 


23 


0 o 2 2 2 2i 

{(n 2 - n 3 + k g2 - k e3 ) + [ 2 ( n 2 k - c3 - n 3 k c2 )3 } 

(n 2 + n 3 ) 2 Uk c2 + k c3 ) 2 


(9) 


2 *1 

‘ 12 " C(ni + n 2 ) 2 + k c 2 2 W 2 * 


( 10 ) 


‘23 


2{[n 2 (n 2 +n 3 ) + k c 2^ k c2 +]< c3^ + ^ n 2 k c3~ n 3 k c2' ) \ 

(n 2 +.n 3 ) + < k c2 +k c 3 ) 


1/2 


( 11 ) 


x _ 2ir , 

6 ■ X V * 


( 12 ) 


* = X k c 2 d * 


(13) 


♦ = Tan’ 1 (— f~ 2' * 


2 k c2 n l 


n l * “2 ‘ k o2 



(14) 
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k 23 


= Tan " 1 { 2(n2k ° 3 


n 3 k c2 ) 


2 2 ‘ . 5 

"2 - "3 - k 03 


+ k 


c 2 


(15) 


The quantities r^j r 23 5 "^ 12 * an< ^ ^23 Can viewed as the magnitudes of 
the complex Fresnal Coefficients. The phase shifts <t >^ 0 and - ^23 are ^ ue 
to absorption in the antireflection layer and the silicon respectively • 

In the above formulations "n" is the real part of the index of refraction 
and k c is the imaginary part. For a transparent . film , k^ is zero, as is 

In performing this calculation with a silicon substrate, experimental 
data for a, n^, and kg versus wavelength were utilized. Results yielding 
the transmission coefficient as a function of photon energy are shown in 
Fig. 5. The reflectance for a bare silicon surface (d = 0) shown in the 
, figure agrees directly with the data of H. R. Phillip [10], Also shown is 
the transmission coefficient for an antireflection film of SiC^ at its 
optimum thickness. For SiO , the transmission, reflection, and 
absorption coefficients are illustrated In Fig. 6 . These optimum thicknesses 
are determined by the maximum efficiency of photon transmission into the 
substrate as described below. 

The results for photon transmission efficiencies for both SiO and 

SiO ^ are shown in Fig. 7 as functions of antireflection layer thickness. 

o 

This illustrates an optimum thickness of around 11G0A for SiOg an< * 

O 

800A for SiO. It has been assumed in this calculation that the SiO 2 
layer is transparent within the spectral region of interest with a 
constant index of refraction of 1.44., This is a good assumption fqr 
SiOgi but SIO is somewhat lossey in this region. To include this 
absorption, the data of Hass and Salzburg [16] of the Index of refraction 


\ 



and extinction coefficient was used to determine values for n ^ and 

k _ for SiO. This data is illustrated in Fig. 8. 
c2 

The available experimental data on reflectance from silicon surfaces 
is somewhat scattered. Prince has reported a maximum reflectance of 6 % 
in the spectrum below 3 ev [4], while Phillip [10] and Malovetskaya fl7] 
present data illustrating reflectances ranging over 30% within the same 
spectral energies. The reflectances calculated in this report for 
bare silicon (Fig. 5) agree closely with the most recent experimental 
data of Phillip. Although the program utilized the index of refraction 
and absorption coefficient data calculated by Phillip from the 
experimental reflectance measurements, it should be realized, however, 
that the absorption coefficient data is essentially the same experimental 
data as that published by Dash and Neumann Cll]« Although Malovetskaya 
does not report . the thickness of the SiO^ layer utilized in his work, 
he does report a maximum transmission of about 93% at 1.5 eV for this 
case. This agrees quite well with the calculations presented in this 


paper . 
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III. CALCULATION TECHNIQUES 

Utilizing the data and equations discussed in the previous 
section^ a computer program has been developed to calculate the carrier 
generation rate under various conditions including that of a general 
single layer antireflection coating. The general flow of the program 
is illustrated in Fig. 9. The incident photon density is determined 
from published data on full spectrum solar irradiance . A most 
comprehensive collection of data on solar radiant energy is the 
NASA document for Space Vehicle Design Criteria which develops final 
data from weighted averages of several investigations [18]. The 
parameters of interest are the solar constant and the solar spectral 
irradiance. The data used for irradiance under AMI conditions (sea 
level) is that of P. Moon [19]. Due to great uncertainties as to 
atmospheric absorption and latitude of measurement this data must be 
viewed as typical. The data was read numerically and stored in the 
program. The incident photon density is related to the photon energy, 
irradiance, and wavelength interval as 

<£> M 2 

N = (#/sec/om ) (16) 

phr 

where (— ) is the average irradiance over an increment of AA , and E 
A A A 

is the average energy of the A A increment. The published data is a 

LJ 

tabular form of ( 77 -) versus wavelength* All available data above 1.11 ^eV 

A A 

and up to about 7 eV was included. 


The next stage in the calculation is to read into the program ■the- 
published data for the absorption coefficient and index of refraction 
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for silicon. This data was supplied in tabulated form from H. R. Phillip 
[10 1 . These values are interpolated to the wavelengths of the spectral 
irradiance data. This then supplies N a s n 3» an ^ ^ or ea °k 
wavelength available from the solar spectrum data. With these values 
the reflectance of the Si surface and the Si-SiC^ surface can be 
calculated. For the SiO layer a similar interpolation process on the 
index of refraction (n^) and extinction coefficient data [16] 

allowed the calculation of these coefficients at the relevant wavelengths . 
At this point the surface photon efficiency is calculated as 


% = 


EH.. T ru 
phx A 
A 


l N 


(17) 


phi 


where T is the transmission coefficient. The surface loss of Table 1 is 


i-v 

Prior to calculating the generation rate, a spatial distribution 
is developed which limits the change in G e between spatial points to 
about 10%. This then allows the generation rates to be calculated at 
each spatial point by summing the contribution to of each incident 
wavelength. Once this data is developed in a tabular form, it is 
numerically integrated to obtain the spatial dependence of the maximum 
optical current . 

The total number of incident photons in the visible spectrum for 

energies above the silicon band edge (1,11 eV for this work) sum to 
17 2 17 2 

3.31 x 10 /cm“ for AMO and 2.14 x 10 /cm for AMI. These values are 
about 5% below what is normally used. There is, however, some ambiguity 
as to the proper value of these quantities since the published irradiance 
data for AMO and AMI vary somewhat. Also, since silicon is an indirect 



band gap material, there is not a sharp edge to the absorption curve making 

the assignment of the band gap energy somewhat ambiguous. For the AMO 

condition a 3% change in the energy band gap (i_.e_. from 1.08 eV to 1.11 eV) 

can result in a 3% decrease in the total number of incident photons. 

Photons at these energies, however, contribute little to the photocurrent 

due to the extremely low absorption coefficients in the immediate vicinity 

of the band edge. The photon density from the irradiance data corresponds 

2 2 

to an available current of 52.9 mA/cm for AMO and 34.2 mA/cm for AMI. 

The inclusion of enchanced quantum efficiency at higher energies resulted in 
a 1,5% increase in the total available current over the total spectral 
current for AMO. The increase for AMI was quite insignificant, being below 



IV. CONCLUSIONS 


It is found that the spatial decay of the optical generation rate 

in silicon due to solar irradiance is extremely rapid. This very rapid 

decay may allow narrower devices in collection applications than 

previously expected; however, it also indicates the importance of 

surface treatments to greatly reduce surface recombination. For AMO 

applications the SiO and Si0 2 antireflection layers are seen to be 

roughly equivalent, with the SiO layer slightly more efficient. This 

rough equivalence points out the need to include the absorption 

characteristics of the SiO film, for this degrades the efficiency of 

this film by roughly 5%. However, the optimum thickness for the SiO 

Q o 

layer is somewhat thinner than the Si0 2> 800A vs . 1100A. In the actual 
use of SiO as an antireflection layer on solar cells the difference 
between SiO and Si0 2 is likely to be greater than indicated by the 
photon efficiency calculation. The SiO film results in a lowered 
generation rate near the surface and an enhanced generation rate deep 
within the silicon. Thus surface recombination losses may be less 
with SiO than with Si0 2 (assuming the same surface recombination 
velocity) giving SiO a slightly larger advantage over Si0 2> For AMI 
applications the optimum film thicknesses are approximately the same 
as AMO. The SiO film is still somwhat more efficient than the Si0 2 
with, the' difference sligh.tly"more : 'pronounc6d' since the AMI spectrum ■■ 
is quite deficient in higher energy photons-. 
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Table 1. Summary of Excess Carrier Generation in Silicon 


Geometry 

v . . 

Spectral 

Conditions 

Optimum Anti- 
Refleotio^ 
Thickness (A) 

Surface 

Loss'* 1 

(%) 

Total Optical 
Current - ** 
(mA/cm 2 ) 

Surface 

Generation 

Rate'*’ 

(#/ec/sec) 

Si 

AMO 

N. A. 

36.4 

34.2 

22 

1.15x10 

AMI 

N. A. 

34.7 

22.4 

21 

1.62x10 

Si + SiO 

AMO 

800 

15.6 

45.4 • 

l 

21 

5.96x10 

AMI 

800 

10.4 

30.7 

21 

1.39x10 

i .... 

Si + SiO ^ 

AMO 

1100 

17.6 

44.3 

22 

1.25x10 

AMI 

1100 

14*5 

29.3 

21 

1.83x10 


■("Computed at Optimum Antireflection Thickness if applicable. 
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Figure 1 
Figure 2 

Figure 3 

Figure 4 
Figure S 
Figure 6 

Figure 7 
Figure 8 


Spatial dependence of the excess carrier generation rate for 
AMI and AM0 illumination. 

Spatial dependence of the excess carrier generation rate for 
both single layer SiO and Si0 2 films at their optimum 
thickness. 

Percentage of optical current collected versus device 
length. Although, this curve is for AMO with no anti- 
reflection coating the' other conditions have quite similar 
results . 

Structure utilized in the calculation of the optical 
coefficients , 

Transmission coefficients for silicon and silicon plus the 
optimum thickness of Si0 2 . 

. Transmission (T) , Absorption (A) , and Reflection (R) 
coefficients for optimum thickness of a single layer SiO 
film. 

. Efficiency of photon transmission into silicon versus 
film thickness for both SIO and SiO^. 

. . Real and imaginary parts of the index of refraction for 
SiO [16]. 


Figure 9. Flow chart for optical generation program. 
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Figure 4. Structure utilised In the calculation of the optical- 
coefficients- 
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Figure a.' Floy: chart for optical generation program. 












Appendix C: Minority Carrier Reflecting Properties of 

Semiconductor High -Low Junctions 


J. R. Hauser and P. M. Dunbar 
Electrical Engineering Department 
North Carolina State University 
Raleigh, North Carolina 27607 

ABSTRACT 

This paper discusses the results of some .computer calculations 
on the minority carrier reflecting properties of high-low junctions. 

In general it has been found that high-low junctions obey first 
order device models to a high degree of accuracy. The ratio of 
carrier densities on the sides of the junction depletion region is 
exponentially related to the high-low junction voltage, which' is constant 
until high injection effects begin to occur. For high injection on the 
lightly doped side of the high-low junction, the junction potential 
decreases and the minority carrier reflecting properties of the 


junction are degraded. 



I . INTRODUCTION 


High-low junctions (either p + -p or n -n) find a number of 
applications in solid state devices. They are most frequently 
used as an interface region between a more lightly doped semiconductor 
and a metal contact in order to provide an ohmic contact to the 
semiconductor . In this application they are known to provide an 
essentially linear current - voltage relationship for majority carrier 
flow across the high-low interface Cl,2j. In addition to aiding in 
making ohmic contacts, high-low junctions also act as minority 
carrier reflecting boundaries . . ' In many minority carrier devices 
this is a useful side effect which can significantly reduce minority 
carrier recombination at ohmic contacts. 

The present work reports on some computer solutions of the 
minority carrier reflecting properties of high-low junctions which 
are located near an injecting p-n junction. This, work has resulted 
from a study of p + -p-n + type diode structures. However, the same 
type of behavior should be observed at any high-low junction 
located near an injecting p-n junction. A comparison is also presented 
of first order analytical models for high- low junctions with the 
exact computer solutions. 
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II. HIGH-LOW JUNCTION PROPERTIES 

The device structure studied in this work is shown m Fig. 1. 

It is an n + -p-p + diode with the n + diffused layer modeled by a 

^ nrt 2Q, 3 

Gaussian impurity profile with a surface concentration of 19 / cm 

and a junction depth of 0*25 pm. The p-region is taken as a uniformly • 

doped region with results for 10 ft* cm and 2000 ft* cm material reported 

in this. work. The p + region has been taken as a uniformly doped 

region with an impurity concentration of 10 /cm and a thickness 

of either 5 pm or 0.5 pm. The overall device thickness is 250 pm 

for the results reported in this paper. 

r The electron lifetime in the p-region was taken as 100 psec and 
1000 psec for the 10 ft*cm and 2000 ft- cm material respectively. Other 
device thicknesses have been studied but the total device thickness 
does not strongly inf luence - the minority carrier reflecting properties 
of the p + -p junction. 

The device structures studied have been less than a diffusion 
length in total thickness. The p + -p interface thus acts as a minority 
carrier reflecting boundary preventing electrons injected by the n + -p 
junction from reaching the back ohmic contact. Various injected electron 
concentration profiles for this diode can appear as shown in Fig. 2 as 
the back boundary ranges from a perfectly reflecting contact to an ohmic 
contact. For the devices studied in this work the high-low junction 
closely approximates an ideal reflecting back contact even for p + 
thicknesses of 0.5 pm. 

There are several parameters which characterize a high-low junction 
as far as its minority carrier reflecting properties are concerned. 



First, there is the equilibrium high-low junction voltage, 
is related to the doping densities as 


, which 


3 


V.. = — CN +/N ) , 

hlo q p p 


Cl) 


where N . and N are the p -1- - and p-region doping densities. The high- 
P 1- P 

low junction potential can only change from its equilibrium value if 

there is a significant change in the majority carrier concentration on 

one or the other side of the high-low junction. A second fundamental 

parameter is the ratio of minority carriers Celectrons m this case) 

across the junction. Assuming that the change in minority carrier 

quasi-Fermi level is negligible across the interface gives Cl J , 

n + = n exp(-qV^/kT) . (2> 

P P 

If high injection does not occur on the lightly doped p-side, V^-is 

approximately the equilibrium value and 

n + a n N /N + . (3) 

P P P P 

This expression is the basis for the reflecting properties of the high- 
low junction since it shows that the minority carrier concentration is 
reduced by the ratio N /N^ + in going from the p to p + side of the 
junction. 

The expression of Equation C3) is based upon a low injection 
condition. In the more general case, Equation (2) must be combined witl) 
the assumption of charge neutrality at the edges of the junction depletion 
region to give the following more complete expression: 

n N 

= n (l + 77^0 — . 

p N- N . 

P P + 


n + 
P 


( 4 ) 
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The additional correction factor of Cl + n /N ) becomes important when 
high injection occurs on the lightly doped side of the high-low junction* 
A final parameter characterizing the minority carrier reflecting 
properties of a high-low junction is what Gunn [3j terms the junction 
leakage velocity or which can also be termed the effective surface 
recombination velocity, S + , high-low junction. As opposed to 

the ratio of minority carriers across the interface which depends only 
on the high-low junction doping levels, the leakage velocity depends 
on the width of the p + region and the minority carrier lifetime in t^is 
region. For uniform doping in the p + region and neglecting recombination 
within the high-low junction space charge region, 


S + 
PP^ 


D . W + n + 

_2Ei cota (Ji, , 

np+ np T p 


(5) 


where D . and L + are the electron diffusion constant and diffusion 
np T np T 

length in the p+ region. The above can be. combined with Equation .(4) 
to give 

■ .V 


s + 

pp 


D + N W + n 

t ' H^ 1 coth (■— ■ £ — ) C 1 + . 

L j. Lj. N 


( 6 ) 


np r p T 


np n 


At low injection this is a constant but the value increases linearly 
with at high injection. This causes the minority carrier current 
into the high-low junction to increase linearly with carrier density 


n at low injection and to increase as n at high injection. 

P P 
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III. COMPUTED RESULTS 


The simple models discussed in fhe previous section have been 
found to agree extremely well with, exact computer calculations for 
high-low junctions * The computer program used in the present device 
analysis has been previously discussed Calculations have been 

made for the structure of Fig. 1 for a variety of dimensions and 
doping densities » Calculations are made at various terminal voltages 
which result in increasing levels of carrier injection into the p-region* 
Shown in Fig. 3 are curves for electron concentration plotted from the 
back surface through the high-low junction for a 0.5 pm wide p + region* 
The electron concentration in the p-side ranges from slightly above 
10 9 /om 3 to about 10 16 /cm 3 while the concentration on the p + side ranges 
from about 10 4 /cm 3 to about 10 14 /cm 3 . The values on the curves are 
the values of applied diode voltage corresponding to each injection 
level. For the low values of the junction depletion region extends 
beyond 1.5 pm and the carrier densities are not yet constant. Shown 


in Fig. 4 is the dependence of.n + on n^ for this device (curve A). 
The solid curve is the simple model as given by Equation .(*+) while 


the data points are results of the exact computer calculations. The 

onset of high injection can be seen for the highest injection levels. 

Calculations for a 5 pm wide p + region are essentially identical 

with the results shown in Fig. 4 with the ratio of minority carriers 

across the high-low junction independent of the width of the p + region. 

Also shown in Fig. 4 are similar data (curve B) for a diode with a 

15 3 

p-side concentration of 1.3x10 /cm . At high injection hoth curves 
approach the same limiting values which are_independent of p-side 




Under high injection conditions, the potential across the high-low 
junction decreases from its equilibrium value. This is shown in Fig. 5 
where the decrease in potential is seen for the four highest injection 
levels. For curves (a) and (b) the junction depletion region extends 
beyond 1.5 pm. As the junction potential decreases at high injection, 
the ratio of minority carriers across the interface decreases and the 
reflecting properties of the high- low junction are degraded. The 
effective surface recombination velocity of the junction increases as 


high injection occurs according to Equation (6). 

Shown in Fig. 6 is a comparison of the high-low junction leakage 
current as given by the expression 


J" , - — qn S j. , 
nhl ^ p pp + 


( 7 ) 


and the exact results of the computer solutions. Again it is seen that 

at low injection the leakage current, is proportional to n^ indicating 

a constant S + while at high injection the leakage current depends on 
PP T 

n^ indicating that s pp + approaches a linear dependence on n^. As • 
indicated in the figure, at low injection the high-low junction has a 
low effective surface recombination velocity. The nearly 
ideal minority carrier reflecting properties of these high- low junctions 
is due to this low value of effective surface recombination velocity; 
Similar high-low junctions but with a 5 pm thick p layer have about an 
order of magnitude smaller effective surface recombination velocity. 
Implicit in the derivation of Equations (3) and (4) is the 


assumption that the minority carrier quasi-Fermi level does not change 
in going across the high-low junction. This assumption has been verified 
with the computer calculations. Shown in Fig, 7 is the calculated 
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quasi-Fermi levels for electrons ^ around the high— low junction* The 

only change in <j>_^ around the high-low junction is the almost abrupt 

change in slope of the <j> curves at the high-low junction. This is 

dcfi n 

to be -expected since ~~ must be inversely proportional to the electron 
concentration in a region where there is negligible change in current 
density. 




a 


IV. CONCLUSIONS 

This work has reviewed "the first order approximations normally 
employed to describe the minority carrier reflecting properties of 
high'- low junctions and compared these expressions with exact computed 
results. The first order results are found to be in very good agree- 
ment with the computer results when account is taken of high injection 
effects on the lightly doped side of the high-low junction. High 
injection reduces the effectiveness of the high-low junction as a 
minority carrier reflecting boundary and leads to an increased 
effective surface recombination velocity of the junction. 
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FIGURE CAPTIONS 

. Device structure for n + -p-p + diode (A) 1000Q* cm p-region, 

N = 7xl0 12 /om 3 ,., N + = 10 1B /cm 3 (B) 10£2»- cm p-region, 

= 1.3xlQ^ 3 /cm 3 ,^N + = lO^/cm 3 . 

P P + 

.. Effect of back contact on injected electron concentration 

(W «L ); (a) Ideal reflecting back contact (b) partially 

p n 

reflecting back contact (c) ohmic contact. 

. Calculated minority carrier concentration around the high -low 
junction at various injection levels (1000 £2* cm p-region). 

, Minority carrier concentration at edge of depletion region 
in p + ~side as a function of minority carrier concentration 
at edge of depletion region on p-side. The solid curves 
are simple calculations from Eq. (4)^ while the data points 
are results of computer analysis, 

. Electrostatic potential around high-low junction at various 
injection levels (1000 £2* cm p-region). 

. High-low junction leakage current as a function of. injection 
level on p-side of junction. The solid curves are simple 
calculations from Eqs. (6) and (7) while the data points are 
results of computer analysis. 

. Electron quasi-Fermi level variation around high-low junction 
(1000 £2* cm p-region). 



Device structure for n -p-p diode (A) 10000* cm p-region 

TO 3 ID 3 

N = 7x10 /cm , N + = 10 /cm (B) 100* cm p-region, 

xt p • , _ ..13, 3 P . T '18/ 3 

N “ 1*3x10 /cm , N + = 10 /cm ♦ 






Fig. 2. 


Effect of back contact on injected electron concentration 
(W <<L ); Ca) Ideal . ref lecting back contact (b) partially 
P n 


reflecting back contact 






1000 rL*cm p region 
10 n. cm p region 


•slope a 2 


/ 7 

i A— slope = I 


Minority carrier concentration at edge of depletion region 
In p + -side as a function of minority carrier concentration 
at edge of depletion region on p-si.de . The solid curves 
are simple calculations from Eq» (4) > while the data points 
are results of computer analysis. 
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000 rx-cm p region 
10 xv cm p region 
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High-low junction leakage current as a function of injection 
level on p-side of junction. The solid curves are.«,simple 
calculations ‘from Eqs, (6) and (7) while the data points 
are results of computer analysis. 
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Fig. 7. 


Electron quasi-Fermi 
(1000 n* cm p-region). 


level variation around high-low junction 




